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elementary  predictions  of  the  magnitude  of  the  effects  can  be  sometimes  modi¬ 
fied  by  other  features  of  the  source  mechanism,  the  geological  variations  along 
the  wave  paths  and  within  the  fault  zone.  Third,  the  ratio  of  peak  horizontal 
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velocities  and  least  for  peak  accelerations.  High  frequency  ground  accelera¬ 
tions  show  variations  due  to  scattering,  attenuation,  and  source  asperities 
that  mask  directivity  effects.  , 

The  conclusions  are  based  largely  on  analyses  of  strong  ground  motions 
recorded  in  the  1979  Coyote  Lake,  1979  Imperial  Valley,  and  1980  Livermore 
earthquakes,  all  in  California.  Each  earthquake  provides  positive  evidence 
for  a  measurable  but  variable  directivity  focusing. 
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PREFACE 


This  report  was  prepared  by  Dr.  Bruce  A.  Bolt,  Seismographic  Station, 
Department  of  Geology  and  Geophysics,  University  of  California,  Berkeley, 
under  contract  DACW39-82-M-1125.  It  is  part  of  ongoing  work  at  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  in  the  Civil  Works 
Investigation  Study,  "Methodologies  for  Selecting  Design  Earthquakes," 
sponsored  by  the  Office,  Chief  of  Engineers,  U.  S.  Army.  The  technical 
monitor  is  Mr.  Paul  R.  Fisher  of  the  Office,  Chief  of  Engineers,  U.  S.  Army 

A  number  of  sections  summarize  and  otherwise  depend  heavily  on  work 
of  Dr.  J.  P.  Singh.  As  indicated  in  the  references,  the  detailed  analysis 
is  available  in  "The  Influence  of  Seismic  Source  Directivity  on  Strong 
Ground  Motions,"  Ph.D.  Dissertation,  U.  C.  Berkeley,  1981.  Special  thanks 
for  discussions  and  other  assistance  are  made  to  Dr.  R.  Uhrhammer,  Mr. 

N.  Abrahamson,  and  Mr.  S.  Blakeslee.  Assistance  with  tables  and  figures 
was  provided  by  Mr.  R.  McKenzie. 

Preparation  of  this  report  was  under  the  supervision  of  Dr.  E.  L. 
Krinitzsky,  Engineering  Geology  and  Rock  Mechanics  Division  (EGRMD) , 
Geotechnical  Laboratory  (GL) ,  and  the  general  supervision  of  Dr.  D.  C. 
Banks,  Chief,  EGRMD,  and  Dr.  W.  F.  Marcuson  III,  Chief,  GL. 

Commander  and  Director  of  WES  during  the  preparation  of  this  report 
was  Colonel  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  Fred  R. 
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PART  I 


THE  PROBLEM 


1.1  Objectives  of  Study 

Major  modern  extensions  of  the  model  of  elastic  rebound  along  a 
rupturing  fault  proposed  by  H.  F.  Reid  following  the  1906  earthquake 
relate  to  the  progression  of  the  dislocation  along  the  fault.  The  move¬ 
ment  of  the  source  immediately  introduces  effects  on  the  generated  seismic 
waves  because  of  the  relative  motion  between  the  dislocation  and  the 
surrounding  rock.  Although  actually  contained  in  Reid's  discussion,  it 
is  only  recently  that  the  problem  of  seismic  wave  generation  by  a  moving 
I  source  has  been  studied  quantitatively,  either  theoretically  or  from  the 

point  of  view  of  the  effects  on  ground  motions  of  interest  to  engineering 
and  structural  design. 

The  main  object  of  this  paper  is  to  describe  available  evidence  for 
effects  of  moving  sources  in  earthquakes  and  to  analyze  them  in  terms  of 
the  available  physics  of  wave  emission.  Based  on  this  survey  and  analysis, 
some  suggestions  are  made  on  the  significance  and  incorporation  of  the 
effect  in  strong  ground  motion  estimation.  The  study  is  addressed  mainly 
to  questions  arising  in  geotechnical  studies  for  engineering  purposes.  The 
importance  of  the  moving  source  effect  has  been  raised  but  little  discussion 
is  available  on  it  in  textbooks;  it  is  hoped  that  this  general  report  may 
be  helpful. 

As  the  theory  of  the  emission  of  sound  from  a  moving  acoustic  source 
indicates, simple  directivity  factors  can  be  derived  in  mathematical  form 
which  take  account  of  the  case  in  which  waves  are  produced  by  a  finite  line 
source  travelling  at  uniform  speed  in  a  fixed  plane.  Even  in  this  case. 


however,  the  waves  produced  are  not,  in  general,  plane  waves  and  the 
mathematical  analysis  is  not  elementary.  In  practical  problems,  where 
comparisons  between  theoretical  models  and  actual  observed  strong  ground 
motions  are  to  be  made,  additional  complications  are  likely  because  of 
variations  of  speed  and  intensity  (e.g.,  changes  in  stress  drop)  along 
the  dislocation  front.  In  actual  earthquakes,  the  fault  zones  are  seldom 
simple  planes  separating  perfectly  elastic  media.  Another  complication 
is  that,  at  least  in  the  case  of  large  magnitude  earthquakes,  the  moving 
source  of  the  earthquake  progresses  within  a  definite  fracture  zone 
produced  by  previous  faulting.  The  new  dislocation  or  re-direction  in 
this  case  occurs  in  a  zone  which  may  resemble  a  seam  of  lower  velocity  and 
higher  attenuating  material  than  the  surrounding  country  rock.  In  research 
work  that  has  grown  out  of  this  study,  source  models  with  such  seams  are 
being  investigated  to  determine  the  effect  of  these  more  realistic  geo¬ 
logical  conditions  on  the  generation  of  the  elastic  wave  field  by  the  moving 
source. 

Finally,  it  is  necessary  to  define  the  words  "focusing"  and  "directivity" 
used  in  the  title  of  this  report.  Wave  focusing  normally  means  the  con¬ 
centration  of  wave  energy  by  variations  in  the  refractive  indices  of  the 
transmitting  media,  as  in  focusing  of  light  by  an  optical  lens  or  focusing 
of  P  and  S  waves  by  changes  in  the  elastic  moduli  at  structural  boundaries. 

The  usage  here  is  different.  It  is  meant  to  describe  an  azimuthal  change 
in  wave  energy  around  the  seismic  source  due  to  the  moving  rupture.  This 
azimuthal  dependence  is  a  function  of  the  direction  of  rupture  and  hence  is 
described  as  a  directivity  effect.  In  theoretical  seismology,  the  ratio 
of  the  spectral  displacements  of  the  seismic  waves  at  two  diametrically 
opposite  points  of  the  fault  plane  is  termed  the  directivity  (see  Part  III). 
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The  accepted  physical  source  model  for  tectonic  earthquakes  is  a 
rupture  which  extends  over  a  fault  plane  in  the  rocks  of  the  Earth  (Fig¬ 
ure  1.1).  The  rupture  progresses  as  a  series  of  dislocations  after 
initiating  at  some  point  (called  "the  seismic  focus").  A  dislocation 
travels  with  various  rupture  speeds  which  are  governed  by  patches  of 
roughness  or  asperities  along  the  fault  itself.  At  the  dislocation  front, 
the  rock  slips  in  a  finite  time,  called  "the  rise  time,"  and  thus  produces 
an  elastic  rebound  of  each  side  of  the  fault.  The  elastic  rebound  reduces 
the  strain  and  hence  the  stress  along  the  fault,  giving  rise  to  local 
"stress  drops."  As  the  released  elastic  strain  energy  is  converted  to 
local  heat  energy  and  wave  energy,  the  seismic  waves  are  generated  near 
the  dislocation  front  and  are  radiated  away  from  the  moving  seismic  source. 

Small  earthquakes  can  be  modeled  as  essentially  point  sources  because 
the  dislocation  front  moves  only  a  few  kilometers  away  from  the  focus. 

On  the  other  hand,  earthquakes  with  engineering  significance  involve  rup¬ 
ture  lengths  of  tens,  and  even  hundreds,  of  kilometers  in  the  greatest 
seismic  sources.  The  appropriate  mechanical  model  must  therefore  take 
account  of  the  movement  of  the  source  over  a  finite  area  with  the  consequent 
superposition  of  radiated  waves  of  various  kinds  from  all  parts  of  the  rup¬ 
tured  front.  In  this  model,  the  duration  of  strong  shaking  is  dependent 
upon  the  dimensions  of  the  faulted  surface.  The  length  of  time  it  takes 
for  a  train  of  seismic  waves  to  reach  a  certain  point  will  depend  upon  the 
total  duration  to  complete  the  dislocation  process  (involving  the  rupture 
velocity)  as  well  as  the  average  seismic  moment  of  the  earthquake  and  the 
size  of  the  stress  drops  along  the  fault. 


RECEIVER 


A  Unidirectionally  Propagating  Rectangular  Strike-Slip  Fault  Model 
Represented  by  Moving  Point  Sources  (After  Singh,  1981) 


Deviations  from  the  mean  duration  from  earthquake  to  earthquake  are 
observed  and  this  is  often  due  to  the  multiplicity  of  the  earthquake 
source.  Many  earthquakes  are  observed  to  have  extended  durations  which 
are  explicable  in  terms  of  the  addition  of  a  number  of  earthquake  ruptures 
in  succession  along  the  fault.  An  important  example  is  the  1940  Imperial 
Valley  earthquake;  a  comparison  of  accelerations  at  the  same  site  between 
the  colocated  1940  and  1979  earthquakes  is  given  in  Figure  1.2.  An  explanation 
is  decrease  in  the  speed  of  the  dislocation  (sometimes  to  zero),  followed 
by  a  readjustment  of  the  elastic  strain  and  the  reinitiation  or  accelera¬ 
tion  of  the  rupture.  The  process  may  be  repeated  so  that  the  full  earth¬ 
quake  is  produced  by  multiple  later  ruptures  from  which  body  waves  propa¬ 
gate  to  the  site  before  the  slower  surface  waves  from  the  previous  ruptures 
have  terminated.  If  the  fault  surface  has  a  large  distribution  of  rough¬ 
ness,  with  asperity  barriers  from  place  to  place  along  the  fault,  then 
considerable  variations  in  source  conditions  might  be  expected  to  arise; 
the  overall  effect  would  be  a  train  of  superimposed  multiple  shocks  in 
rapid  succession.  This  feature  of  many  large  earthquakes  complicates  the 
recorded  strong  motions  and  makes  it  difficult  to  isolate  the  effects,  if 
any,  of  single  aspects  of  the  rupture  process  such  as  directivity  focusing. 


IMPERIAL  VAllEY  EARTHQUAKE 

amax=0.40  g  NORTH  COMPONENT 


A  Comparison  of  Strong  Motion  Records  Obtained  at  the  El  Centro 
Station  (North  Component)  during  May  18,  1940,  and  October  15 
1979,  Earthquakes 


1.3  Available  Evidence 

Strong  azimuthal  patterns  in  seismic  wave  radiation  have  been  observed 
instrumentally  in  both  local  earthquakes  and  teleseisms  for  many  years. 

From  seismic  source  theory,  these  radiation  patterns  have  been  defined  for 
all  types  of  seismic  waves:  P  and  S  waves,  and  Love  and  Rayleigh  waves. 

These  radiation  patterns  are  the  basis  of  the  method  of  fault  plane  solu¬ 
tion  using  the  polarity  (or  phase)  of  first  motions  developed  by  P.  Byerly 
in  the  1930's.  The  Byerly  model,  however,  considers  the  earthquake  source 
to  be  stationary.  When  the  moving  dislocation  is  taken  into  account,  major 
changes  can  occur  in  the  radiation  patterns  for  the  various  seismic  waves. 

A  considerable  amount  of  observational  evidence  consistent  with  such 
modifications  has  now  been  published,  although  most  evidence  pertains  to  motions 
recorded  a  considerable  distance  from  the  seismic  source,  i.e.,  in  the  far 
field. 

The  first  strong  evidence  for  the  effect  of  a  moving  source  in  seis¬ 
mology  was  discussed  by  H.  Benioff  (1955)  in  his  explanation  of  the  in¬ 
tensity  pattern  observed  in  the  1952  Kern  County  (California)  earthquake. 

The  argument  was  that  a  result  of  the  propagation  of  the  dislocation  along 
the  fault  would  be  different  signals  received  at  opposite  ends  of  the 
fault,  with  larger  intensities  of  higher  frequency  in  the  direction  of 
propagation  and  smaller  intensities  and  lower  frequencies  at  the  opposite 
end  (see  Figure  1.3).  Early  seismological  theoretical  developments  were 
made  by  Ben-Menahem  (1961),  who  developed  a  function  for  the  directivity 
based  on  the  known  theoretical  results  for  a  moving  source  in  acoustics 
(see  Appendix  B).  Ben-Menahem  noted  that  the  spherical  shape  of  the  Earth 
provides  a  simple  way  to  apply  directivity  functions  to  the  determination 
of  source  parameters  such  as  fault  length  and  average  rupture  velocity. 
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DIRECTION  OF  PROPAGATION 


Figure  1.3  Illustration  of  Sound  Waves  Proceeding 
from  a  Moving  Source 

Ben-Menahem  examined  the  multiple  arrivals  of  long-period  Love  and  Rayleigh 
waves,  which  had  circled  the  Earth  a  number  of  times  in  the  case  of  several 
major  earthquakes.  In  particular,  he  showed  that  directivity  effect  was 
present  in  the  great  Chilean  earthquake  of  22  May  1960  and  the  Alaskan 
earthquake  of  28  March  1964.  A  recent  example  of  this  type  of  comparison 
in  the  far  field  between  theory  and  observation  will  be  given  in  Section 
4.2. 

There  is  also  available  more  recent  seismological  evidence,  both  from 
the  analysis  of  isoseismal  lines  relative  to  the  ruptured  fault  and  from 
analyses  of  accelerograms  recorded  near  to  small  and  moderate  earthquakes 
in  California  with  well-determined  seismic  sources  (Bakun  et  al.,  1978). 

It  must  be  stressed,  however,  that  much  of  the  evidence  in  both  categories 
is  of  a  circumstantial  nature;  the  recordings  can  be  explained  in  part  by 
directivity  focusing,  but  the  explanation  is  not  unique  in  most  cases. 


Observed  seismic  intensity  distributions  are  seldom  very  precise  and 
are  much  affected  by  variations  in  soil  and  other  surficial  conditions. 

A  review  of  some  of  the  field  evidence  has  been  given  by  J.  P.  Singh  (1981) 
The  theoretical  analyses  with  speculative  but  plausible  models  depend 
significantly  on  assumptions  made  concerning  both  the  seismic  source 
mechanism  and  the  crustal  structure. 

In  geotechnical  work  which  endeavors  to  predict  ground  motion  at  a 
site  by  the  specification  of  strong  ground  motion  parameters,  questions 
on  directivity  focusing  have  been  raised  seriously  only  in  the  last  few 
years.  For  such  purposes,  convincing  evidence  must  be  available  from 
strong-motion  accelerometers  in  the  near  field,  rather  than  from  sensitive 
seismographs  remote  from  the  source  that  record  small  ground  motions. 

At  the  present,  there  are  only  a  few  such  cases  (and  these  from  California) 
where  instrumental  evidence  is  available  on  the  question  of  whether  sig¬ 
nificant  directivity  effects  occur  close  to  the  seismic  source.  The 
most  valuable  cases,  discussed  in  Part  IV,  Involve  the  1979  Coyote  Creek, 
the  1979  Imperial  Valley,  and  the  1980  Livermore  Valley  earthquakes. 


Consequences 


The  problem  of  the  significance  of  directivity  focusing  in  the 
prediction  of  strong  ground  motion  as  a  basis  for  engineering  design 
has  arisen  already  in  practice.  The  thrust  of  the  argument  is  that, 
for  near  fault  sites,  unless  effects  of  directivity  are  included  in 
ground  motion  estimation,  the  results  may  under-estimate  the  actual  maxi¬ 
mum  amplitudes  of  the  ground  shaking.  An  important  illustration  which 
will  be  referred  to  at  several  places  in  this  report  relates  to  arguments 
put  before  the  Nuclear  Regulatory  Commission  that  directivity  be  included 
in  seismic  analysis  for  the  Diablo  Canyon  Nuclear  Power  Plant  site,  on 
the  coast  of  central  California.  This  site  is  only  4-5  miles  from  the 
Hosgri  fault  on  which  a  Richter  magnitude  (Ml)  7.5  earthquake  was  postulated. 

It  should  be  emphasized  that  seismic  directivity  focusing  is  of  geo¬ 
technical  interest  only  as  a  special  near-f ield  earthquake  ground  motion 
effect  which  may  be  experienced  at  a  site  towards  which  fracture  propa¬ 
gation  progresses.  In  theory,  the  effect  of  such  focusing  is  to  cause 
local  amplification  of  the  ground  motion  associated  with  the  earthquake. 

It  is  suggested  that  this  may  be  a  crucial  consideration  in  structural 
design  because  such  focusing  might  significantly  increase  the  spectral 
components  of  ground  motion  employed  as  input  to  the  seismic  design  of 
critical  structures.  In  the  case  of  the  site  of  the  Diablo  Canyon  nuclear 
plant,  it  was  argued  that  a  major  earthquake  (Ml  ■  7.5)  associated  with 
rupture  along  the  Hosgri  fault  zone  may  produce  unexpectedly  high  accel¬ 
erations  in  certain  directions. 

In  its  discussion  of  the  question,  the  Safety  and  Licensing  Appeal 
Board  (Document  No.  50-2750L  and  50-3230L)  even  raised  quite  specific 


issues  and,  for  example,  stated:  "Interveners  and  the  applicant  have 
suggested  that  the  strong-motion  data  obtained  from  stations  along  the 
direction  of  the  Imperial  fault  (in  the  1979  Imperial  Valley  earthquake) 
evidence  the  'focusing  of  earthquake  motion.'  Yet  when  the  acceleration 
data  of  two  such  stations,  El  Centro  Arrays  Nos.  6  and  7,  are  plotted  as 
a  function  of  distance  from  the  fault,  the  horizontal  acceleration  values 
fall  well  below  the  regression  mean  line  for  1  km  distance.  The  vertical 
acceleration  values  are  also  lower  than  the  mean  on  such  a  fault.  To  the 
extent  possible,  the  parties  should  analyze  the  seismic  records  of  this 
earthquake  as  they  pertain  to  the  focusing  phenomena  and  relate  the 
results  of  such  analyses  to  the  likelihood  that  focusing  might  result  in 
amplified  seismic  motion." 

It  is  now  generally  accepted,  based  on  the  earthquake  source  models 
discussed  above,  that  directivity  focusing  occurs  in  every  earthquake  to 
some  degree.  It  is  still  a  matter  of  debate,  however,  whether  in  most 
cases  the  seismic  source  is  organized  and  coherent  enough  to  affect 
materially,  by  means  of  source  directivity,  all  high-frequency  amplitudes 
(e-g-t  peak  ground  accelerations)  as  a  function  of  azimuth  (particularly 
concentrated  in  a  narrow  cone  along  the  ruptured  fault) . 

Further  practical  aspects  of  the  matter  related  to  assessment  of 
ground  motions  for  engineering  purposes  and  hazard  analysis  will  be  dis¬ 
cussed  in  more  detail  in  Section  5.2  after  the  various  lines  of  evidence 


are  reviewed. 


According  to  some  theoretical  estimates  of  the  effects  of  directivity 
focusing,  the  amplitudes  of  seismic  waves  travelling  ahead  of  the  moving 
dislocation  might  be  greater  by  an  order  of  magnitude  than  those  propa¬ 
gating  in  the  opposite  direction  of  the  dislocation.  Such  theoretical 
results  must  be  interpreted  in  terms  of  the  frequencies  of  the  waves  and 
the  rupture  conditions.  Because,  for  example,  all  ruptures  propagate  up¬ 
wards  from  the  focus,  to  a  certain  extent  every  major  earthquake  contains 
the  effect  of  directivity  focusing  at  points  on  the  surface  where  the 
rupture  breaks  out. 

Another  important  aspect  is  the  occurrence  of  simultaneous  dislocation 
in  various  directions  along  the  fault  surface.  It  is  known,  for  example, 
that  the  fault  source  in  many  earthquakes  ruptures  in  a  bilateral  manner, 
with  the  rupture  extending  away  from  the  focal  region  in  opposite  di¬ 
rections.  An  example  of  bilateral  rupture  is  the  San  Andreas  dislocation 
in  the  1906  San  Francisco  earthquake  where  rupture  commenced  somewhere  off 
the  Golden  Gate  of  San  Francisco  Bay  and  propagated  both  to  the  north  and 
to  the  south  for  hundreds  of  kilometers.  In  many  other  earthquakes,  it 
has  been  established  that  rupture  is  predominantly  unidirectional;  examples 
may  be  found  in  damaging  Turkish  earthquakes  along  the  Anatolian  fault 
(Richter,  1938)  and  in  the  1966  Parkfield  earthquake. 

With  these  source  complexities  in  mind,  intensity  studies  from  great 
earthquakes  provide  general  information  on  extent  and  significance  of  any 
amplification  due  to  directivity  focusing.  In  the  case  of  the  1906  San 


Francisco  earthquake,  where  detailed  field  reports  are  available  on  damage 
including  many  photographs,  there  is  little  evidence  for  enhanced  damage 
due  to  source  directivity.  Although  the  isoseismal  map  published  in  the 
Report  of  the  State  Earthquake  Investigation  Commission  shows  high  inten¬ 
sities  along  the  San  Andreas  fault  break,  critical  reassessments  by  G.  D. 
Louderback  and  others  demonstrate  that  damage  to  structures  was  not  par¬ 
ticularly  great  along  the  fault  zone  in  this  earthquake  (Bolt,  1981). 
Louderback  not  only  pointed  to  this  fact  but  also  suggested  that  the  reason 
may  well  be  the  presence  of  gouge  and  shattered  rocks  in  the  fault  zone 
(see  Section  3. A). 

Another  example  of  a  great  earthquake  which  involved  surface  rupture 
of  a  well-marked  fault  is  the  4  February  1976  Guatemala  earthquake  produced 
by  rupture  along  the  Motagua  fault  (see  Figure  2.1).  Espinosa  (1976)  has 
drawn  attention  to  the  asymmetry  in  the  seismic  intensity  as  shown  by  the 
isoseismal  lines  drawn  in  Figure  2.1.  Rupture  progressed  from  the  focus 
near  Los  Amates  at  the  east  end  of  the  fault  to  west  of  Guatemala  City. 

A  much  larger  area  of  Modified  Mercalli  (MM)  intensities  VIII  and  VI  was  re¬ 
ported  at  the  western  portion  of  the  meizoseismal  region  and  this  is  consis¬ 
tent  in  a  general  way  with  the  concept  of  directivity  focusing.  It  should  be 
noted  that  isolated  pockets  of  intensity  IX  and  VIII  were  mapped  at 
various  places  along  the  fault.  Field  investigations  show  that  at  places 
along  the  fault  rupture  damage  to  certain  types  of  ordinary  buildings  was 
not  great. 

Finally,  in  any  study  of  seismic  intensity,  in  making  any  inferences 
based  on  field  observations,  it  must  be  emphasized  that  damage  to  struc¬ 
tures  is  likely  to  arise  predominantly  from  high*frequency  motions  (in 
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Although  many  Isoselsmal  maps  for  earthquakes  in  the  United  States 
and  elsewhere  have  been  published,  few  are  based  on  field  observations 
with  the  precision  necessary  to  make  a  decisive  statement  on  the  effective¬ 
ness  of  directivity  focusing.  Below  a  few  of  the  most  definitive  are 
mentioned. 

One  example  which  led  to  seismological  speculation  on  the  importance 
of  source  focusing  was  the  earthquake  of  21  July  1952  (Ml  =  7.2)  in  the 
series  of  earthquakes  in  Kern  County,  California.  In  this  shock,  fault 
rupture  began  at  a  focus  about  16  km  below  the  surface  along  the  White  Wolf 
fault.  It  is  believed  that  the  rupture  extended  upward  and  westward, 
breaking  the  surface  along  the  surface  trace  (see  Bolt,  1981).  The  in¬ 
tensity  maps  do  not  plot  detailed  isoseismals  in  the  meizoseismal  area 
around  the  fault  itself  where  intensities  were  VIII  or  over.  Richter  (1958) 
mentions  the  possible  confusion  in  assignment  of  intensity  arising  from  the 
rupture  effects  at  the  fault  trace  itself,  a  similar  circumstance  to 
Louderback's  point  of  the  1906  earthquake.  Richter  points  out,  "There 
was  no  evidence  of  extremely  violent  shaking  at  points  near  the  trace,  as 
judged  by  houses.  Intensity  IX  was  manifested  over  much  of  the  area  near 
the  fault  and  effects  due  to  shaking  assignable  to  X  were  developed  in  many 
localities."  Another  feature  of  the  isoselsmal  map  in  this  case  was  the 
general  elongation  of  the  outer  isoseismals  in  the  northwest -southeast 
direction  which  is  parallel  to  the  regional  structure,  although  the  White 
Wolf  fault  itself  trends  almost  at  right  angles  to  this  structure. 

On  the  other  hand,  in  this  earthquake  at  the  longer  periods  there  was 
directivity  focusing  at  teleseismic  distances.  As  quoted  in  Richter's  book, 
"It  was  found  that  the  amplitude  of  surface  waves  for  European  recordings 
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were  unexpectedly  high.  However,  the  amplitudes  recorded  in  Australia 
and  New  Zealand  were  correspondingly  low;  at  those  stations  the  second 
group  of  surface  waves  which  had  traveled  around  the  world  by  the  major 
arc  was  actually  larger  than  the  direct  waves  across  the  Pacific.  These 
results  could  be  explained  in  terms  of  piling  of  seismic  wave  energy  in 
the  direction  in  which  rupture  progressed  from  the  instrumentally  deter¬ 
mined  hypocenter  along  the  White  Wolf  fault;  the  speed  of  progression  of 
such  rupture  might  be  expected  to  be  of  the  same  order  as  that  of  surface 
waves.  Something  related  to  the  Doppler  effect  would  result,  with  increased 
concentration  of  energy  in  the  direction  of  faulting;  this  might  have 
contributed  to  the  apparent  intensity  at  Tehachapi."  As  pointed  out  at 
the  end  of  Section  2.1,  however,  the  last  comment  does  not  necessarily 
follow  from  the  directivity  evidence  at  distant  stations,  because  at  the 
teleseismic  distances,  surface  waves  have  longer  periods  than  the  ground 
motions  that  damage  small  structures  such  as  those  at  Tehachapi. 

Other  examples  of  moderate  earthquakes  in  which  isoseismals  are 
available  are  the  1966  Parkfleld  and  the  1971  San  Fernando  earthquakes 
(see  Bolt,  1981).  In  these  cases,  Brune  (1978)  has  stated,  "Directivity 
was  probably  very  important  in  generating  the  high  velocities  observed  for 
the  Parkfleld  and  San  Fernando  earthquakes.  ...  The  concept  of  focusing  or 
directivity  is  important  in  strong  motion  seismology  because,  in  addition 
to  the  effect  of  radiation  patterns,  it  can  introduce  large  azimuthal 
variation  in  ground  motion  that  introduces  a  large  range  of  scatter  in  the 
data,  thus  making  it  particularly  difficult  to  estimate  the  mean  and  stand¬ 
ard  deviation  of  the  expected  velocities  and  accelerations  from  a  limited 
set  of  data  close  to  the  source." 


Summary  comments  on  the  instrumental  evidence  in  the  case  of  these 
two  earthquakes  are  given  in  Section  4.1.  As  far  as  the  influence  of 
directivity  on  the  isoseismals  is  concerned,  the  evidence,  however,  is 
not  so  clearcut.  The  1971  San  Fernando  earthquake  had  a  combined  thrust 
and  left-lateral  fault  source.  The  intensity  distribution  for  this  earth¬ 
quake  is  shown  in  Figure  2.2.  The  highest  intensity  area  (VII)  extends 
from  the  fault  source  towards  the  southeast  in  an  elongated  oval  pattern. 
It  is  aligned  with  the  direction  of  slip,  rather  than  parallel  to  the 
strike  of  the  fault,  in  general  agreement  with  the  notion  that  greater 
energy  was  propagated  in  the  direction  of  the  moving  dislocation. 

Two  other  intensity  maps  are  reproduced  here  (Figures  2.3,  2.4)  for 
the  mainshocks  of  the  sequences  near  Coyote  Lake  (6  August  1979)  and 
Livermore  (24  January  1980).  Although  instrumental  measurements,  dis¬ 
cussed  in  detail  in  Sections  4.1  and  4.3,  indicate  significant  source 
focusing  effects  in  the  near  field  of  these  earthquakes,  it  will  be  seen 
that  the  isoseismal  patterns  do  not  provide  much  definitive  evidence  one 
way  or  the  other.  (Information  on  geological  structure  and  soils  would 
also  have  to  be  incorporated  in  any  analysis.)  There  is,  however,  a 
slight  indication  of  extensions  of  the  isoseismals  to  the  southeast  of 
the  foci  in  both  cases. 
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Figure  2.3  isoseismal  map  for  the  central  California  earthquake  of  6  August  1979 
17  05  22.7  UTC.  Roman  numerals  represent  Modified  Mercalli  inten¬ 
sities  between  isoseismals;  Arabic  numerals  are  used  to  represent 
these  intensities  at  specific  sites. 
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Figure  2.4  Isoseismal  map  for  the  central  California  earthquake  of  24  January 
1980,  19  00  09.7  UTC.  Roman  numerals  represent  Modified  Mercalli 
intensities  between  isoseismals;  Arabic  numerals  are  used  to  rep¬ 
resent  these  intensities  at  specific  sites. 


PART  III 


THEORETICAL  BASIS  FOR  SOURCE  DIRECTIVITY 

3.1  Physical  Models 

The  full  detail  of  the  patterns  of  seismic  waves  cannot  be  predicted 
in  the  near  field  unless  the  earthquake  source  is  modeled  in  a  realistic 
way.  The  first  specification  of  the  physics  of  the  generation  of  such 
waves  was  made  by  H.  F.  Reid  in  his  studies  of  the  faulting  that  occurred 
along  the  San  Andreas  fault  in  the  1906  San  Francisco  earthquake.  Reid 
stated  that  the  earthquake  vibrations  "originate  in  the  surface  of  fracture 
the  surface  from  which  they  start  has  at  first  a  very  small  area,  which 
may  quickly  become  very  large,  but  at  a  rate  not  greater  than  the  velocity 
of  compressional  elastic  waves  in  the  rock."  This  extract  from  Reid's 
description  of  the  elastic  rebound  theory  postulates  the  propagation  of 
rupture  along  the  fault.  Further,  Reid  stated,  "It  is  probable  that  the 
whole  movement  at  any  point  did  not  take  place  at  once,  but  proceeded  in 
regular  steps."  It  is  of  interest,  therefore,  that  despite  this  foresight 
definite  consequences  of  the  moving  seismic  source  were  not  observed  instru- 
mentally  in  earthquakes  until  seismograms  of  seismic  surface  waves  at 
distant  stations  were  compared  after  the  great  1960  earthquake  in  Chile. 

In  the  mid-1960s,  N.  Haskell  made  key  contributions  to  the  development 
of  a  more  physically  realistic  model  and  he  added  details  to  the  way  that 
waves  were  generated  along  the  fault.  He  constructed  a  model  "in  which  the 
fault  displacement  is  represented  by  a  coherent  wave  only  over  segments  of 
the  fault  and  the  radiation  from  adjacent  segments  is  assumed  to  be  stat¬ 
istically  independent  or  incoherent."  The  physical  situation  is  that  the 
rupture  begins  suddenly  and  then  spreads  with  periods  of  acceleration  and 


retardation  along  the  relatively  weakly  welded  fault  zone.  A  simplified 
model  of  this  circumstance  is  shown  in  Figure  3.1. 

A  few  years  later  this  Reid-Haskell  model  was  further  refined  by  the 
introduction  of  mechanical  concepts  such  as  areas  of  roughness  or  "as¬ 
perities"  along  the  fault  rupture  surface.  This  representation  has  been 
developed  by  K.  Aki  and  others  (see  Aki  and  Richards,  1980)  who  have 
further  drawn  attention  to  the  possible  importance  of  seismic  "blocks"  or 
barriers  along  the  fault  at  which  higher  strains  are  required  for  the 
moving  dislocation  to  propagate. 

It  must  be  emphasized  that  there  is  difficulty  in  making  the  model 
of  the  last  paragraph  geologically  realistic.  The  estimation  of  strong 
ground  motion  at  any  point  depends  not  only  on  the  ability  to  specify  a 
realistic  source,  but  also  on  the  ability  to  specify  a  realistic  crustal 
structure  in  which  the  rock  layers  have  anelastic  and  structural  properties 
(including  damping)  appropriate  for  the  region  in  question.  Most  theoretical 
work  to  date  assumes  parallel  horizontal  layers  with  soils  underlain  by 
rocks  of  increasing  rigidity.  Sloping  layers  with  the  presence  of  low- 
velocity  zones  and  anomalous  rock  lenses  may  significantly  modify  the  sur¬ 
face  ground  motions  through  well-known  wave  behavior,  such  as  geometrical 
focusing,  diffraction  and  scattering  of  the  seismic  waves.  In  the  1979 
Imperial  Valley  earthquake  (see  Section  4.1),  high-frequency  waves  (par¬ 
ticularly  on  the  vertical-component  records  at  many  sites)  have  been  inter¬ 
preted  as  P  waves  refracted  sharply  upwards  through  surficial  sedimentary 
layers  with  strong  velocity  contrasts.  These  properties  of  the  wave  propa¬ 
gation  path  can  be  dominant  and  make  the  interpretation  of  a  strong-motion 
accelerogram  difficult  and  non-unique.  Unfortunately,  detailed  crustal 
structure  is  not  known  for  most  sites,  with  consequent  additional  un  cer- 


Simplified  Model  of  Rupture  Surface  during  February  9,  1971  San  Fernando 
Earthquake  in  Relation  to  the  San  Fernando  Surface  Faulting,  Pacoima  Dam 
and  the  Focus  (modified  from  Bolt  1972). 


tainties  in  the  prediction  of  a  strong  ground  motion. 

An  Aki-Haskell  earthquake  source  with  a  high  degree  of  physical 
realism  would  be  the  preferred  starting  point  for  the  computation  of 
synthetic  seismograms  of  ground  motion  near  to  the  source  so  that  direct¬ 
ivity  focusing  could  be  correlated  between  prediction  and  observation.  A 
difficulty,  in  practice,  is  that,  in  nearly  all  cases  the  detailed  asperity 
and  barrier  structure  along  the  fault  is  not  known  ab  initio.  Only  after 
the  earthquake,  can  (non-unique)  models  be  contructed  to  fit  the  records. 

Another  critical  matter  is  not  emphasized  particularly  in  this  report. 
The  physics  of  earthquake  sources  indicates  that  wave  patterns  generated 
by  moving  fault  ruptures  might  be  significantly  different  in  certain 
aspects,  depending  on  the  different  types  of  fault  rupture  involved. 

Purely  strike-slip  motion  is  likely  to  produce  in  practice  rather  different 
effects  on  directivity  focusing  than  dip-slip  dislocations. 


3.2  Kinematics 


A  stationary  seismic  fault-slip  source  can  be  modeled  as  a  double¬ 
couple  at  a  point.  It  follows  that  there  is  already  an  azimuthal  radiation 
pattern  in  wave  emplitudes  even  before  a  rupture  velocity  is  introduced. 

An  example  of  this  azimuthal  variation  is  shown  in  Figure  3.2.  For  a 
simple  right-lateral  strike-slip  fault,  for  example,  the  radiation  pattern 
in  the  near  field  for  S  waves  will  involve  four  equal  lobes  with  maxima 
occurring  in  the  direction  of  the  fault  and  at  right  angles  to  it.  This  azi 
muthal  effect  is  not  related  to  a  source  propagation  and  will  always  be 
present  for  each  type  of  seismic  wave  considered.  Movement  of  the  source 
produces  deformation  of  the  lobes  with  asymmetrical  effects  in  front  of 
and  behind  the  moving  dislocation.  Brune  (1978)  states,  "The  concept  of 
focusing  or  directivity  is  important  in  strong  motion  seismology  because,  in 
addition  to  the  effect  of  radiation  pattern,  it  can  introduce  a  large  azi¬ 
muthal  variation  in  ground  motion." 

For  a  moving  source,  waves  that  leave  the  traveling  rupture  in  opposite 
directions  will  generally  have  different  amplitudes.  The  kinematic  or 
geometrical  effect  is  simply  demonstrated  by  means  of  circular  wavefronts 
which  are  emitted  at  equal  intervals  of  time  from  the  moving  point.  The 
case  of  a  moving  point  traveling  along  a  straight  line  is  illustrated  in 
Figure  3.3  which  was  discussed  initially  in  a  seismologicai  context  by 
Benioff  (1955)  in  connection  with  the  1952  Kern  Country  earthquakes. 

The  effect  of  the  movement  of  the  source  is  to  change  the  spectral 
amplitudes  of  the  waves  which  leave  the  source  in  opposite  directions.  The 
ratio  of  the  spectral  amplitudes  is  called  the  directivity  or  directivity 
function  (see  Appendix  B).  The  formula  for  the  directivity  is  easily 
derived  from  simple  kinematic  arguments  and  is  given  in  equation  (3.1): 
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Figure  3.2  Radiation  patterns  of  the  P-  and  S-waves  from 

sources  with  rupture  propagation  for: 

(a)  v  /  v  =0.5;  (b)  v  /v  =0.9,  where  v  and 
it  s  r  s  r 

v  denote  the  rupture  and  shear  wave  velocities 

respectively.  Solid  lines  indicate  positive, 

and  dashed  lines  negative  values. 

(From  Kasahara,  1981) 
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where  b  is  fault  length,  c  is  wave  velocity,  v  is  rupture  velocity 
and  0q  is  azimuthal  angle. 

The  formula  assumes  that  unilateral  faulting  occurs,  that  is  to 
say, the  fracture  along  the  fault  propagates  in  only  one  direction.  If 
the  fracture  moves  in  both  directions,  as  in  bilateral  faulting 
(see  Section  3.1),  then  the  directivity  function  is  more  complicated. 


For  the  sake  of  illustration,  consider  a  vertical  strike-slip 
fault  with  unilateral  fracture  propagation.  The  computed  radiation 
patterns  of  P  and  S  waves  are  shown  in  Figure  3.2.  Here  the  wave 
amplitudes  are  plotted  in  the  x-y  plane  in  the  case  of  the  unilateral 
rupture  moving  in  the  positive  x  direction.  Rupture  velocity  ratios 
are  assumed  in  case  (a)  v/c  =0.5  and  in  case  (b)  v/c  =  0.9  for  the  shear 
wave  velocity  in  the  medium.  (The  latter  ratio  is  more  realistic.)  The 
effect  of  the  moving  source  is  immediately  evident,  since  for  a  stationary 
point  source  the  lobes  would  show  a  symmetrical  radiation  pattern. 

As  discussed  in  the  book  by  Kasahara  (1981),  the  energy  of  the 
seismic  wave  pulse  is  approximately  constant  for  different  azimuths  so 
that  in  the  positive  direction  of  the  x  axis  where  the  amplitude  becomes 
very  large,  the  pulse  width  becomes  quite  narrow. 

It  should  be  emphasized  that  the  discussion  given  above  is  a 
geometrical  one  and  assumes  simple  geometric  configurations.  In  terms 
of  propagation  theory,  the  waves  are  assumed  to  have  phase  coherence 
throughout  their  path,  an  assumption  which  may  be  approximately  so  if 
smooth  rupture  takes  place.  In  many  earthquakes,  however,  as  already 
mentioned,  thete  will  be  "chattering"  of  the  rupture  through  patches 
of  roughness  on  the  fault  plane  and  the  phase  coherence  will  become 
very  degraded. 


3.3  Results  from  Acoustics 

Consider  a  source  of  sound  moving  along  a  path  (Figure  3.4)  specified 

by  the  position  vector  r  (t)  with  components  x  (t),  y  (t)  and  z  (t). 

s  s  s  s 

The  observation  point,  0,  is  at  a  position  r  with  coordinates  x,  y,  z 

at  rest  with  respect  to  the  medium.  The  sound  pressure  observed  at 

r  at  time  t  was  emitted  by  the  source  at  time  t  =  t  -  (R/C)  at  which  time 

the  source  is  at  rest  r  (t  )  at  the  emission  point  E.  The  distance 

s  e 

R  =  r  -  rg(te)  between  the  emission  point,  E,  and  the  observation 
point,  0,  is  determined  by 

R2  =  [x-xg(t-R/ C) ] 2  +  [y-ys(t-R/C>]2  +  [z-zs(t-R/C)]2  (3.2) 


More  than  one  value  of  R  satisfies  equation  (3.2),  depending  on 
the  Mach  number,  i.e.,  the  ratio  of  the  speed  of  source,  V,  to  the 
speed  of  sound,  C.  There  are  two  cases:  subsonic  or  supersonic  which 
are  considered  separately. 

Subsonic  case  (M  ■  V/C  <1) 

In  the  simplest  case  of  subsonic  motion  the  source  moves  along 

the  x-axis  (see  Figure  3.4).  if  the  source  passes  the  origin  of  the 

coordinate  system  at  time  t  =  0,  then  x  (t)  =  Vt  and  Y  =  Z  =0.  Then 

s  s  s 

equation  (3.2)  simplifies  to 

R2  =  [x-V  (t-R/ C)  ] 2  +  r2  =  [(x-Vt+MR]2  +  r2  (3.3) 

2  2  2 
where  r  *  y  +  z  . 

Equation  (3.3)  is  satisfied  by 


For  M<1,  only  the  positive  sign  gives  positive  R. 


Figure  3.4  Kinematics  of  a  point  source  in  subsonic  motion. 


For  M>1,  both  plus  and  minus  in  equation  (3.4)  are  allowable.  Positive 


R  is  taken  when  the  point  of  observation  is  behind  the  source,  i.e., 
(x-Vt<C).  Then, 

_  _  M(Vt-x)  ±  (Vt-x) J l-(M2-l)r2/ (Vt-x)^ 


Now,  let  r/ (Vt-x)  =  tan  0^  and  sin  0^  =  1/M  as  shown  in  Figure  3.5, 
yielding 

M(Vt-x)  ±  (Vt-x)  /l-(tan0  /tan0  )2 

R+  - - 5 - - - —  •  (3.6) 

M  -1 


R+  is  real  only  when  ©2<0l  so  Point  of  observation  is  contained 

inside  a  cone  with  the  vertex  at  the  source  and  the  vertex  angle  20^. 

In  this  case,  two  emission  points  E+  and  E_  can  be  constructed  with  emission 
times  corresponding  to  R+  and  R_  that  produce  simultaneous  contributions 
to  the  sound  wave  field  at  the  point  of  observation  0  at  time  t.  Thus 
point  0  lies  on  the  intersection  between  two  phase  surfaces  (Figure  3.5) 
which  are  both  tangents  to  the  cone  and  have  their  centers  (E+  and  E_) 
on  the  x-axis. 

For  the  degenerate  case  where  the  observation  point,  0,  is  located 
along  the  x-axis,  the  distinction  between  the  kinematics  for  subsonic 
and  supersonic  sources  can  be  visualized  in  a  simple  manner  in  a  space- 
time  diagram  (Figure  3.6).  In  the  subsonic  case  a  point,  0,  will  be 
reached  at  a  certain  time  only  by  the  forward -going  wave  from  the  source 
corresponding  to  the  emission  point  E.  In  the  supersonic  case,  there  is 
one  forward-  and  one  backward-going  wave  reaching  0  simultaneously, 
although  emitted  at  different  times  corresponding  to  E  and  E_. 


Supersonic,  V>C 


Subsonic,  V<C 


Figure  3.6  Space»time  diagrams  for  the  wave  kinematics  of 

sources  in  uniform  subsonic  and  supersonic  motions. 

The  kinematic  properties  for  a  moving  acoustical  source  in  terms 
of  frequency  and  amplitude  modulation  may  now  be  developed.  First, 
consider  a  monopole  source  of  sound  moving  with  uniform  velocity  with 
respect  to  the  surrounding  medium.  The  monopole  acoustic  source  can  be 
thought  of  in  terms  of  a  pulsating  sphere,  the  strength  of  which  is 
specified  by  a  scalar  q. 

For  the  subsonic  case,  the  azimuthal  variation  of  the  sound  pressure 
from  this  moving  monopole  source  can  be  expressed  (Morse  and  Ingard,  1968) 
as 

p  =  _L  <i'  (t-R/C)  +  q (cos  0-  M)V 

4ttR  (1-M  cosO)^  4ttr2(1-M  cos©)^  '  (3.7) 

where  the  prime  denotes  differentiation  with  time  and  the  first  term 
represents  the  radiation  field  with  corresponding  pressure  decrease  as 
R  ^ .  Note  that  when  0=90°,  the  pressure  in  the  far  field  is  the  same  as 
when  the  source  is  at  rest.  Quantitative  analysis  of  equation  (3.7)  shows 
that  the  contribution  to  the  total  pressure  field  from  the  second  term  in 
the  equation  is  small  (generally,  less  than  about  2-3%).  Therefore, 
dropping  the  second  term,  the  approximate  pressure  field  is 


p  = 


1 

4ttR 


q'  f(t— (R/C)] 
(1-M  cosG)" 


(3.8) 


This  relationship  shows  that  the  sound  pressure  in  the  direction  of  motion 
(0=0°)  is  larger  than  the  pressure  in  the  direction  away  from  the  direction 
of  motion  (©=180°)  by  a  factor  (1+M) (1-M) ^ . 

Now  let  the  source  strength  be  harmonic,  i.e.,  q(t)  =  qQ  sin  u >Q(t). 

The  term  q'  inequation  (3.9)  becomes  qoL0Qcos  [a)Q (t-R/C) ] .  The  pressure 
field  then  has  the  phase 


<j>  =  oi  (t-R/C) 
o 


Note  that  this  phase  <}>  is  not  proportional  to  time  alone  because  R  also 
is  time  dependent.  A  generalized  concept  of  frequency  can,  however,  be 
defined  as  the  time  derivative  of  the  phase,  i.e. 


0)  = 


d<{> 

dt 


„  1  <1]U 

“o  (1"  C  dt5 


Substitute 


1  dR 
c  dt 


-  M  cos  0 
1-M  cos  0 


yielding 


0)  =  (D  (1+ 

o 


M  cos  0  . 
1-M  cosQ' 


1 


GO 


o 

-M  cos  0 


(3.9) 


Equation  (3 . 9)  is  the  Doppler  formula  for  frequency.  As  the  moving 
source  passes  the  observing  point,  the  angle  0  increases  from  0=0  to 
0=ir  so  that  the  frequency  increases  from  goq/(1-M)  to  w  /  (1+M) . 

The  Doppler  effect  is  related  to  the  observed  intensity  distribution 
of  sound.  A  point  in  front  of  the  source  receives  the  energy  emitted 
during  the  time  At  in  the  time  interval  At^  which  is  smaller  than  the  time 


interval  At-  corresponding  to  a  point  of  observation  behind  the  source. 


This  result  demonstrates  that  a  source  which  has  a  symmetric  radiation 

pattern  at  rest  focuses  radiation  in  the  direction  of  motion.  If  A  ,  w 

o  o 

and  Aq  are  the  wavelength,  frequency  and  amplitude  for  the  source  at  rest, 
then  for  the  case  of  a  moving  source  the  wavelength  frequency  and  amplitude 
will  be  modified  as  shown  in  Table  3.1  and  Figure  3.7. 


TABLE  3.1 

MODIFICATIONS  TO  WAVE  PROPERTIES 


Forward  Direction 

Backward  Direction 

(0  =  0°) 

(0  =  180°) 

<D 

0) 

Frequency 

a),  =  -q — —  (increase) 
t  1-M 

ub  =  Tim  (decrease) 

A 

A 

Amplitude 

A  =  - -  (increase) 

A  =  - —  (decrease) 

(1-M)Z 

“  (14Mr 

Wavelength 

Af  =  \q  (1-M)  (decrease) 

A.  =  A  (1+M)  (increase 
b  o 

Duration 

Decrease 

Increase 

Because  for  seismic  sources  such  as  fault  dislocations  the  observed  rupture 
velocities  appear  to  be  less  than  the  shear  wave  velocity  in  the  strained 
rocks  along  the  fault ,  the  supersonic  rupture  case  is  not  pursued  further 
here  (but  see  Section  5.3). 

One  interesting  feature  of  the  supersonic  case,  however,  is  the 
interference  among  wave  contributions  from  the  two  emission  points 
E+  and  E  .  These  waves  with  different  frequencies  and  amplitudes  give 
rise  to  amplitudes  of  sound  pressure  that  oscillate  with  time. 

Because  of  the  mathematical  equivalencies,  in  appropriate  cases, the 
above  theory  for  a  moving  acoustical  source  can  be  related  to  P  and  S 
waves  and  surface  waves  from  a  moving  seismic  source.  In  seismology  also, 
the  frequencies  and  the  amplitudes  of  the  seismic  wave  field  will  be 
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modulated  by  geometrical  factors  which  are  a  function  of  a  corresponding 
seismic  Mach  number  M  and  the  azimuthal  angle  6  between  the  earthquake 


source  and  recording  site.  The  usual  factor  sin  x  in  the  spectrum  of 
an  observed  strong  motion  pulse  shape  is  equivalent  in  the  time  domain 
to  a  convolution  with  a  rectangular  function  of  temporal  duration 

cos  9, 


b  (1/V  - 


This  gives  the  apparent  duration  of  rupture  detected  by  seismographs 
along  the  direction  0q  for  rupture  speed  V  and  wave  propagation  speed  C. 
At  higher  frequencies,  the  envelope  of  x  ^  sin  x  is  proportional  to  w  \ 
This  smoothing  is  weakest  in  the  direction  of  rupture  propagation  (0^  =  O' 
and  strongest  in  the  opposite  direction  (0^  =  tt).  As  a  result,  there  are 
more  high-frequency  waves  at  0^  =  0  than  at  0^  =  ir.  This  frequency  var¬ 
iation  is  similar  to  the  Doppler  effect,  which  shifts  the  frequency  u>  of 
a  moving  oscillator  by  [1  -  V/C  (cos  0^)]  In  the  case  of  *  1  sin  x, 
the  nodes  are  shifted  as  a  function  of  0^  as 

a)  .  =^22^  (1-(V/C)  cos  0  )  n-1,  2,  ... 

nodes  b  o 


The  result  is  destructive  interference  between  waves  being  generated  at 
different  parts  of  the  fault  plane;  this  tends  to  smooth  out  high  fre¬ 
quencies.  In  contrast,  the  Doppler  effect  does  not  smooth,  because  there 
Is  only  one  oscillator. 

It  should  be  remarked  that  in  addition  to  the  patterns  produced  by 
the  factor  x  ^  sin  x  there  may  be  other  significant  factors  such  as  a 


superposition  of  sources  from  multiple  complex  ruptures  and  effects  from 
travel  path  variations.  Similar  problems  of  interference  also  arise  for 
moving  acoustical  sources.  Morse  and  Ingard  (1968)  note  that  near  to  the 


3.4  Seismological  Results.  Seam  Waves 


As  was  pointed  out  in  the  description  of  seismic  source  models  given 
in  Section  3.1,  most  numerical  modeling  of  the  earthquake  source  process 
has  assumed  a  moving  dislocation  on  a  single  plane  meant  to  represent 
the  plane  of  the  slipped  fault  in  the  crust.  Geological  conditions, 
however,  along  major  faults  like  the  San  Andreas  are  significantly  dif¬ 
ferent  with  usually  a  rather  broad  zone  (up  to  500  meters  or  more  in  hori¬ 
zontal  width)  of  crushed  and  sheared  rock  often  containing  bands  of  gouge 
extending  to  considerable  depth.  The  question  arises  as  to  what  extent  theo¬ 
retical  predictions  for  models  with  a  single  fault  plane  separating 
perfectly  elastic  media  are  relevant  to  the  prediction  of  ground  motion 
in  a  large  earthquake  very  close  to  the  causative  fault  zone.  Examples  °f 
ite  zones  occur  with  the  San  Andreas  fault  in  California  and  its  main  subsidi¬ 
aries  such  as  the  Calaveras  Hayward  fault  in  northern  California  and  the  San 
Jacinto  fault,  the  White  Wolf  fault,  and  the  Imperial  Valley  fault  in  southern 
California.  To  the  present,  little  numerical  modeling  is  available  to 
test  the  influence  of  such  a  zone  on  the  seismic  waves  generated  from  a 
dislocation  moving  within  it. 

There  would  be  two  effects  of  the  zone.  First  is  the  effect  to  the  fault 

zone  structure  itself.  Geometrically,  it  resembles  a  slab  or  "seam"  like  a 

coal  seam  through  country  rock.  For  the  sake  of  clarity  in  describing  this 

structure,  the  term  "fault  seam"  will  be  used  herein.  The  defining 

properties  of  the  seam  would  be  lower  seismic  velocity  for  the  rocks 

within  the  seam  relative  to  those  outside  the  seam.  There  would  also  be 
a  higher  attenuation  with  the  seismic  parameter  Q  being  up  to  20  percent 
lower  on  the  average  than  the  Q  values  in  the  unfractured  rock  outside  the 


I 

seam.  In  practice,  of  course,  gradations  must  be  expected  between  the 
j  properties  of  the  seam  and  the  rock  more  remote  from  the  major  faults 

in  the  seam  and  also  variations  in  properties  along  the  fault  zone  within 
the  seam.  For  the  purpose  of  exploration,  however,  these  more  complex 
|  effects  will  not  be  taken  into  account  at  the  present  time. 

As  the  seismic  waves  are  generated  by  the  moving  dislocation  at 
depth,  it  is  expected  that  some  of  the  energy  would  be  channeled  within 
j  the  seam  and  that  these  channel  waves  would  be  refracted  out  of  the  seam; 

seismic  wave  energy  would  leak  into  the  country  rock,  dependent  on  the 
lengths  and  types  of  the  waves  in  question.  The  effect  would  happen  to 
j  body  waves  of  both  P  and  S  type  and  to  both  Rayleigh  and  Love  waves  or 

their  equivalents.  There  would  be  partitioning  in  the  wave  energy  between 
the  various  types  of  waves  as  they  interact  with  the  seam  boundaries,  and 
special  types  of  trapped  waves  with  appropriate  wave  lengths  might  become 
established  (analogous  to  Stoneley  waves  which  are  known  to  travel  under 
certain  circumstances  along  boundaries  between  rocks  of  different  prop- 
|  erties  at  depth  in  the  Earth.)  Of  special  importance  would  be  enhanced 

attenuation  of  waves  which  were  trapped  within  the  seam;  wave  energy  dif¬ 
fracted  out  of  the  seam  channel  would  attenuate  less  from  damping  than 
wave  energy  within  the  attenuating  seam  itself.  This  effect  might  be 
expected  particularly  to  decrease  the  high-frequency  accelerations  associ¬ 
ated  with  velocities  and  displacement  (see  Singh,  1981). 

Theoretical  work  supporting  these  ideas  may  be  found  in  the  physics 
literature,  particularly  that  dealing  with  the  channeling  of  light, 
acoustic  energy,  and  elastic  waves  along  ridges  and  ducts  (Morse  and  Ingard, 
1968).  As  a  simple  model,  there  is  some  resemblance  to  the  well-known 
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effect  of  an  aperture  on  a  light  beam  in  which  diffraction  patterns  are 
set  up  (see  Appendix  A,  Figure  A.l).  Singh  (1981)  has  drawn  attention 
to  the  correspondence  between  the  amplitude  pattern  on  the  far  side  of 
a  light  aperture  and  that  which  he  observed  in  the  variation  of  the  peak 
ground  acceleration  on  the  array  across  the  Imperial  fault  in  the  1979 
earthquake  (see  Figure  4.7). 

Some  numerical  analogies  come  from  the  study  of  wave  propagation 
along  irregular  coal  seams  (Krey,  1963;  Drake  and  Asten,  1982).  It  is 
known  that  channel  or  seam  waves  (see  Appendix  A,  Figure  A.l)  propagate 
in  the  two-dimensional  wave  guide  which  is  constituted  by  a  seam  of  low- 
velocity  coal  between  higher  velocity  rock  layers.  These  seam  waves  are 
both  longitudinal  (i.e.,  sound)  waves  and  shear  waves  which  can  be  studied 
by  seismological  methods.  Theory  indicates  that  a  low-velocity  layer 
supports  two  classes  of  seam  waves,  generalized  Rayleigh  waves  and  gen¬ 
eralized  Love  waves.  In  the  latter,  the  motion  is  horizontal  shear, 
transverse  to  the  direction  of  propagation  along  the  seam.  Each  class  of 
seam  waves  includes  multiple  modes  which  are  dispersive. 

Recent  detailed  finite  element  modeling  of  waves  channeling  along 
coal  seams  have  been  carried  out  by  Drake  and  Asten.  They  show  that  waves 
essentially  of  Stoneley  type  can  channel  along  coal  seams  in  both  pseudo- 
Rayleigh  and  pseudo-Love  types  but  that  the  efficiency  of  the  channel  ef¬ 
fect  is  dependent  significantly  on  the  wave  mode  and  on  the  frequency  of 
the  wave.  Typical  P  velocities  in  the  country  rock  in  this  modeling  were 
taken  as  4  km/sec  and  2.1  km/sec,  respectively,  while  within  the  seam  the 
P  and  S  velocities  were  2.4  and  1.28  km/sec,  respectively.  Although  these 


models  do  not  Include  any  damping  terms  in  the  finite  element  computation, 
they  do  show  that  the  presence  of  a  fault  zone  with  significantly  different 
seismic  parameters  could  have  significant  consequences  as  far  as  direct¬ 
ivity  focusing  properties  are  concerned. 

Another  relevant  study  involves  the  propagation  of  elastic  surface 
waves  guided  by  ridges  on  the  surface  of  an  elastic  body.  The  wave  guiding 
action  of  the  linear  topographic  feature  is  relevant  to  the  present  problem 
in  which  we  consider  the  fault  zone  or  fault  seam  as  the  wave  guide  analogy 
of  the  raised  ridge.  Some  useful  results  (Burridge  and  Sabina,  1972)  for 
such  models  using  finite  element  methods  show  that  with  certain  geometric 
dimensions  and  elastic  constants,  guided  waves  do  propagate  along  the  ridge 
as  dispersed  trains.  The  calculations  indicate,  however,  that  in  any  case 
only  some  modes  may  be  properly  trapped  by  the  linear  feature  and  that  this 
occurs  in  only  certain  limited  frequency  ranges.  Again,  the  calculations 
indicate  that  only  some  modes  of  wave  propagation  of  pseudo-Rayleigh  type 
propagate  in  an  unattenuated  form.  In  other  words,  energy  is  continually 
leaked  from  the  wave  guide  into  the  surrounding  material  in  most  cases. 

By  analogy,  for  the  fault  seam  the  amplitude  of  the  waves  propagating  along 
the  seam  would  continually  fall.  As  with  the  coal  seam  analysis,  this 
interesting  modeling  work  was  carried  out  without  any  frictional  attenuation 
included  in  the  calculations,  so  the  results  can  be  taken  only  as  suggestive 
It  is  likely,  however,  that  the  inclusion  of  damping  would  leave  the  main 
results  unaffected  but  reduce  even  more  the  amount  of  high-frequency 
seismic  energy  able  to  be  propagated  to  any  distance  along  the  seam. 
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In  the  last  few  years  Instrumental  evidence  of  Importance  has  come 
from  a  number  of  moderate  California  earthquakes.  These  are  the  1966 
Parkfield  earthquake,  the  1971  San  Fernando  earthquake,  the  1979  Coyote 
Lake  earthquake,  the  1979  Imperial  Valley  earthquake,  and  the  1980  Livermore 
earthquake.  In  each  case,  strong-motion  accelerometers  recorded  the  strong 
shaking  in  the  near  field  and  so  some  quantitative  evidence  is  available 
which  bears  on  the  question  of  the  radiation  wave  field  from  a  moving  dis¬ 
location  (see  for  fuller  descriptions.  Bolt,  1981).  As  is  usually  the 
case,  however,  in  strong-motion  seismology  the  geographical  distribution 
of  strong-motion  instrument  sites  was  by  no  means  optimum  for  the  assess¬ 
ment  of  the  quantitative  effects  of  directivity  focusing.  While  the 
results  are  helpful,  in  no  case  is  there  a  final  conclusion  to  the  problem. 

The  Parkfield  earthquake  occurred  on  28  June  1966  and  had  ML  =  5.5 
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and  a  seismic  moment  of  1.9  x  10  dyne-cm.  The  mean  stress  drop  in  this 
earthquake  has  been  estimated  at  quite  low  values  of  the  order  of  5  bars. 
Rupture  occurred  along  a  central  segment  of  the  San  Andreas  fault  showing 
predominant  strike-slip  motion  with  a  maximum  slip  of  about  20-25  cm. 

The  rupture  was  right-lateral  and  there  is  evidence  of  a  mean  rupture 
velocity  of  between  2.5  and  2.9  km/ sec.  A  small  array  of  strong-motion 
instruments  at  right  angles  at  the  southern  end  of  the  fault  rupture  re¬ 
corded  the  motion  and  has  given  information  on  the  attenuation  of  strong 
shaking  away  from  the  causative  fault.  However,  no  strong-motion  instrument 
was  available  at  the  northern  end  of  the  faulting  so  a  direct  calculation 
of  the  directivity  is  not  possible. 


It  is  of  interest,  however,  that  the  strike  of  the  fault  is  such  that 
its  northward  projection  passes  close  to  the  Berkeley  stations  in  the  San 
Francisco  Bay  area.  A  study  of  the  spectrum  of  the  recorded  waves  on  a 
broadband  seismograph  at  Berkeley  (Filson  and  McEvilly,  1973)  has  given  an 
excellent  example  of  an  observed  directivity  spectral  function  which  fits 
well  the  form  of  the  theoretical  spectrum  for  a  moving  source.  This 
spectrum  yields  an  estimate  of  the  velocity  of  rupture  of  about  2.2  km/ sec 
towards  the  south. 

The  San  Fernando  earthquake  occurred  on  9  February  1971  and  had  a 
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magnitude  of  Ml  =  6.5  and  a  moment  of  Mq  =  1.5  x  10  dyne-cm.  Various 
stress  drops  were  originally  calculated  for  this  earthquake  with  an  average 
value  of  50  bars.  In  a  careful  reassessment,  however.  Hanks  and  McGuire 
(1981)  inferred  that  the  major  part  of  the  seismic  energy  was  released  over 
a  small  portion  of  the  fault  plane  and  that  the  stress  drop  involved  was  of 
the  order  of  1000  bars. 

Surface  faulting  occurred  in  the  San  Fernando  Valley  and  along  the 
foothills  of  the  San  Gabriel  Mountains  in  this  earthquake.  It  indicated 
both  thrusting  and  left-lateral  motion.  One  explanation  given  for  large 
accelerations  recorded  at  Pacoima  Dam  (up  dip  along  the  thrust  fault)  is 
directivity  focusing  (Heaton  and  Helmberger,  1979;  McGuire  and  Hanks,  1980). 
A  number  of  plausible  but  competing  models  have  been  computed  for  this 
earthquake  source,  however,  so  that  available  data  do  not  provide  the 
resolution  needed  for  a  well-constrained  solution. 

The  Coyote  Lake  earthquake  occurred  on  6  August  1979.  It  had  ML  =  5.9 
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and  a  seismic  moment  of  6  x  10  dyne-cm.  Minor  breaks  and  cracks  occurred 
along  the  Calaveras  fault  for  about  8  to  10  km,  with  right-lateral  dis¬ 
placement  of  up  to  about  5  mm.  It  is  of  special  value  for  studies  of 
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directivity  focusing  because  24  strong-motion  accelerometers  were  triggered 
located  within  50  km  of  the  fault  rupture  with  several  very  near  to  the 
Calaveras  fault.  The  geometry  of  the  faulting  and  the  location  of  after¬ 
shocks  and  accelerometers  are  shown  in  Figures  4.1  and  4.2.  Of  particular 
interest  is  a  comparison  (Figure  4.3)  of  recorded  ground  motion  between 
the  Coyote  Creek  station  to  the  north  and  Station  6  to  the  south,  both 
adjacent  to  the  fault.  Comparisons  of  peak  acceleration  (PGA),  peak 
velocity  (PGV),  and  peak  displacement  (PGD)  in  these  records  have  been  made 
in  detail  by  J.P.  Singh  (1981).  He  found  that  enhanced  acceleration, 
velocity,  and  displacement  recorded  at  Station  6  compared  with  the  Coyote 
Creek  station  are  in  agreement  with  directivity  focusing.  The  enhancement 
is  greatest  for  PGD  and  least  for  PGA  in  this  comparison,  with  the  ampli¬ 
fication  factor  reaching  1.8  for  PGA,  2.0  for  PGV,  and  4.0  for  PGD. 
Moreover,  higher  frequencies  were  richer  at  Coyote  Creek  than  at  Station  6 
(Figure  4.4),  perhaps  indicating  damping  in  the  fault  zone. 

The  Imperial  Valley  earthquake  occurred  on  15  October  1979.  It  had 
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a  magnitude  of  Ml  =  6.6,  and  a  seismic  moment  of  8.7  x  10  dyne-cm.  As 
in  the  Parkfield  earthquake,  the  estimated  mean  stress  drop  in  this  earth¬ 
quake  is  low  at  about  5.5  bars.  Rupture  was  observed  (see  Figure  4.5) 
along  the  35  km  of  the  Imperial  and  Brawley  faults  and  the  motion  was 
right-lateral  strike-slip  (up  to  55  cm)  with  some  dip-slip  offsets  (up  to 
19  cm)  down  to  the  east.  Many  strong-motion  accelerometers,  including  a 
13-element  linear  array,  recorded  this  earthquake  in  the  near  field  and 
provided  much  hard  data  on  the  surface  radiation  pattern.  As  mentioned  in 
Section  3.2,  theoretical  modeling  of  wave  focusing  due  to  moving  sources 
suggests  that  the  zone  of  wave  enhancement  is  concentrated  into  a  narrow 
lobe  in  the  direction  of  the  moving  fault  dislocation  (see  Figure  3.2  ). 


Numbered  circles  correspond  to  Table  S-1  and  their 
diameter  scales  to  the  local  magnitude  (M^) . 

Inset  shows  the  relative  position  of  the  Coyote 
Lake  area  within  the  central  coast  region  of 
California. 


Figure  4.1  Coyote  Lake  Earthquake  Sequence  (Modified 
from  Uhrhammer,  1980) 
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Figure  4.5  Location  Map  of  Strong  Motion  Instruments  (modified 
from  Brady  et  al.,  1980) 
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Simple  models  (Brune,  1978)  suggest  that  a  focusing  lobe  has  a  width  of  ±5% 
or  so  from  the  direction  of  the  dislocation.  In  such  a  circumstance,  the 
restricted  lobe  of  enhanced  directivity  focusing  could  be  missed  ob- 
servationally  unless  there  are  a  large  number  of  strong-motion  instruments 
in  the  near  field.  In  the  case  of  the  1979  Imperial  Valley  earthquake, 
with  rupture  traveling  from  near  the  Mexican  border  northwest  beyond 
El  Centro,  the  accelerograms  from  Bonds  Corner  and  Station  7  (see  Figure 
4.6)  are  within  this  lobe  so  that  they  provide  a  test  of  the  predictions. 

The  13-element  array,  total  length  45  km,  passing  through  El  Centro  pro¬ 
vides  further  discriminatory  near-field  evidence. 

The  1979  Imperial  Valley  measurements  have  been  extensively  analyzed 
by  J.P.  Singh  (1981)  and  only  the  main  relevant  conclusions  are  summarized 
here.  At  Bonds  Corner  the  recorded  PGA,  PGV,  and  PGD  (230°  horizontal  com¬ 
ponent)  were  770  cm/sec2,  44  cm/sec,  and  15  cm,  respectively.  At  Station  6 
to  the  northwest,  the  corresponding  values  were  428  cm/ sec2,  108  cm/ sec,  and 
55  cm.  The  latter  values  are  closely  confirmed  by  readings  of  454  cm/sec2, 
108  cm/sec  and  41  cm  at  Station  7.  In  addition,  the  bracketed  duration  at 
Bonds  Corner  (230°  component)  was  13  sec  compared  with  8  sec  at  Station  6. 
This  comparison  indicates  larger  peak  amplitudes  in  velocity  and  displacement 
ground  motions  to  the  northwest  (a  factor  of  about  3)  consistent  with  the 
directivity  model  but  a  decrease  in  PGA  in  the  direction  of  the  rupture 
propagation. 

Another  germane  comparison  is  a  plot  of  PGA,  PGV  and  PGD  values  from 
the  El  Centro  array  (230°  component)  against  perpendicular  distance  from 
the  Imperial  fault  (Figure  4.7).  The  PGV  and  PGD  curves  show  a  roughly 
unimodal  curve  with  maximum  values  extending  1-2  km  to  the  west  of  the 
fault  and  5-6  km  to  the  east.  Outside  this  central  maximum,  wave  attenuation 
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VELOCITY  ACCELERATION  DISPLACEMENT  VELOCITY 


PGV  normalized  to  EL  CENTRO  6 
PCD  normalized  to  EL  CENTRO  6 


PERPENDICULAR  DISTANCE  FROM  IMPERIAL  FAULT  IN  km 


is  high  and  more  or  less  monotonic  reaching  about  one-fifth  maximum 
value  at  a  distance  of  about  20  km.  The  PGA  curve,  by  contrast,  has  a 
broader  central  maximum  and  is  tri-modal.  Stations  2  and  11  at  about 
15  km  from  the  ruptured  fault  show  higher  accelerations  than  expected  from 
the  average  trend.  The  explanation  of  these  "anomalous"  readings  could 
be  unusual  site  soil  conditions  or  unusual  surficial  crustal  properties. 
Singh  (1981)  prefers  to  explain  the  PGA  contours  in  terms  of  a  diffraction 
pattern  set  up  by  a  zone  of  low-velocity  rocks  across  the  fault  zone. 

Seismic  refraction  data  suggest  that  such  a  low-velocity  seam  (see  Section 
3.4)  is  over  1  km  in  width  in  the  basement  rocks.  The  hypothesis  is  that 
the  high-frequency  waves  (5-10  Hz)  recorded  in  the  accelerograms  are  damped 
markedly  along  the  Imperial  fault  zone,  whereas  the  longer  period  waves  in 
the  velocity  and  displacement  records  are  showing  mainly  directivity  effects 
because  the  channeling  and  damping  effects  of  the  fault  seam  are  not  so 
great  at  periods  of  1  sec  and  greater. 

A  final  point  made  by  Singh  relates  to  the  amplitude  spectra  of 
recorded  ground  motions.  Spectral  amplitudes  for  the  230°  horizontal 
component  of  ground  displacement  at  Stations  6,  7  and  Bonds  Corner  are 
plotted  in  Figure  4.8.  The  spectral  envelopes  for  Stations  6  and  7  (to 
NW  of  the  rupture)  exceed  the  spectral  envelope  for  Bonds  Corner  for  fre¬ 
quencies  less  than  0.8  Hz  and  are  less  for  higher  frequencies.  This  result 
is  in  line  with  the  above  argument. 


4.2  Foreign  Earthquakes  \ 

Descriptions  of  seismic  intensity  in  a  number  of  earthquakes  that 
have  been  studied  closely  in  the  field  outside  the  United  States  have  led 
to  suggestions  of  focusing  effects.  Unfortunately,  in  many  cases,  the  use 

of  the  word  "focusing"  is  not  precisely  defined.  The  investigators  sometimes 
have  in  mind  the  focusing  of  waves  by  variations  in  underground  structures 
and  sometimes  the  effects  of  the  dynamics  of  a  moving  source.  One  important 
recent  example  is  the  1977  Romanian  earthquake  (m^  =  6.8).  The  isoseismal 
lines  were  clearly  asymmetric,  extending  to  south  and  southeast  -  a  pattern 
(Bolt,  1981)  noticed  in  earlier  Romanian  earthquakes  from  the  Vrancea  zone 
near  the  Carpathians. 

Another  clear  example  of  a  great  foreign  earthquake  which  produced 
intensities  which  probably  can  best  be  explained  by  directivity  focusing 
is  the  Guatemala  earthquake  of  4  February  1976.  This  earthquake  was  pro¬ 
duced  by  rupture  from  east  to  west  of  the  Motagua  fault  and  the  isoseismal 
map  is  given  in  Figure  2.1.  Modified  Mercalli  Intensities  near  Guatemala 
City  and  to  the  north  were  notably  stronger  than  intensities  in  the  vicinity 
of  Los  Amates,  with  a  spreading  out  of  the  Isoseismal  curves.  Espinosa 
(1976)  is  of  the  opinion  that  some  geometrical  amplification  effects  of 
the  seismic  waves  are  involved. 

As  well,  there  is  direct  seismological  evidence  that  (at  least  at  long 
periods)  directivity  focusing  was  important  in  the  source  mechanism  of  this 
earthquake.  Espinosa  shows  the  amplitudes  of  the  very  long  period  SH  sur¬ 
face  waves  called  G  waves  which  were  recorded  on  long-period  seismographs 
at  the  WWSSN  Station  Porto  in  Portugal.  As  shown  in  Figure  4.9,  the  train 
of  waves,  called  G2» traveling  westward  (i.e.,  in  the  direction  of  the 


source  velocity) around  the  great  circle  from  Guatemala  to  Porto  has  a  large 
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NOTES: 

1)  Waves  C2  and  C4 
left  the  source  in 
direction  of  rupture 
propagation  at  an 
angle  of  about  15°  from 
the  direction  of  fault 
rupture. 

2)  Wave  G3  left  the 
source  at  an  angle  of 
about  165°  from  the 
direction  of  fault 
rupture. 

3)  Wave  C4  has 
traveled  213°  farther 
than  Wave  G3. 


Figure  4.9  Effect  of  Source  Propagation  on  G-Wave  Amplitudes 

Recorded  at  Station  Porto,  Portugal  (from  Espinosa, 

1976) 

hundred-second  wave  pulse.  This  pulse  is  still  evident  in  the  wave  train 
called  G4  which  is  produced  after  the  G2  train  travels  another  circuit 
around  the  Earth.  In  comparison,  the  wave  train  called  G3, which  is  the 
SH  surface  wave  pulse  that  has  left  the  source  in  the  opposite  direction 
to  the  source  propagation, has  a  virtually  undetectable  amplitude  at  Porto. 

Specific  evidence  of  directivity  focusing  has  also  been  given  in  some 
seismological  studies  of  Japanese  earthquakes.  As  yet,  however,  no  cases 
provide  evidence  of  the  same  strength  as  the  California  studies.  It  is 
interesting  that  in  his  recent  book,  K.  Kasahara  (1981),  in  dealing  with 
the  mechanics  of  earthquakes,  discusses  directivity  focusing  as  a  likely 
significant  phenomenon  but  does  not  quote  any  studies  from  Japanese  earth¬ 
quakes  . 
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A  number  of  directivity  functions  have  been  evaluated  in  the  far  field 
for  long-period  Love  and  Rayleigh  waves  and  the  results  published.  These 
results  are  consistent  with  the  effectiveness  of  the  term  given  as  equation 
B-6  in  Appendix  B.  The  examples  include  the  great  Chilean  earthquake  of 
22  May  1960  and  the  great  Alaskan  earthquake  of  28  March  1964  (Ben-Menahem 
and  Singh,  1981). 


4.3  The  Greenville  (Livermore)  Earthquake  Sequence,  24  January  1980 


The  principal  earthquake  at  1900^  00m  09.46s  UTC  on  the  24th  of 
January,  1980,  with  Ml  =  5.5  at  Berkeley,  caused  a  surprising  amount  of 
damage  north  of  Livermore  Valley  and  was  associated  with  surface  rupture 
along  the  Greenville  fault  (Figure  4.10).  There  was  a  foreshock  (Ml  =  2.7) 
a  minute  and  a  half  earlier  and  a  sequence  of  59  events  (ML  _>  2.5)  in  the 
ensuing  six  days.  On  January  27  at  02^  33m  35.96s,  a  second  principal 
earthquake  occurred  in  the  sequence  (Ml  =  5.6  at  Berkeley).  This  second 
principal  shock  was  located  14  km  to  the  south  of  the  first  principal  earth¬ 
quake  towards  the  southern  end  of  the  Greenville  fault.  The  location  of 
the  earthquake  sequence  in  relation  to  the  Greenville  fault  and  the  Livermore 
Valley  in  central  California  is  shown  in  Figure  4.10.  The  two  principal 
earthquakes,  1/24  and  1/27,  are  number  2  and  11,  respectively. 

Preliminary  estimates  of  the  seismic  moments  of  the  two  principal 
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shocks  are  5.3  x  10  and  1.3  x  10  dyne-cm,  respectively.  The  sequence 
is  of  particular  interest  for  the  present  study  for  two  reasons.  First, 
clear  surface  faulting  was  observed  on  a  mapped  fault,  the  Greenville  fault, 
and  extended,  as  observed  by  cracks  in  fields  and  roads,  along  the  strike 
for  at  least  6  km.  The  fault  offset  was  predominately  right-lateral  but 
varied  from  place  to  place  in  the  soil,  sometimes  accompanied  by  comparable 
vertical  offsets  with  small  sag  sections  common.  On  the  average,  the 
northeast  side  went  up.  The  location  of  the  focus  and  aftershocks  argue 
for  rupture  from  north  to  south  along  the  fault  in  the  1/24  mainshock. 

Secondly,  the  radiation  patterns  in  intensity  for  the  two  principal 
shocks  were  clearly  different,  with  the  shock  of  January  27  exhibiting  an 
asymmetrical  intensity  radiation  pattern  with  larger  amplitudes  towards  the 
northeast.  The  first  principal  earthquake  caused  damages  of  various  extents 
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Figure  4.10  Epicenters  of  Principal  Earthquakes  in  the 
Greenville  Sequence,  January  1980 


to  a  few  facilities  of  the  Lawrence  Livermore  National  Laboratory.  Because 
of  its  moderate  size  (Ml  =  5.5),  this  circumstance  caused  some  surprise 
and  raised  questions  of  the  mechanism  of  the  wave  generation  in  this  case. 

The  1/27  principal  earthquake  had  its  focus  near  the  south  end  of  the  Green¬ 
ville  fault.  It  produced  less  intensity  in  the  Livermore  area  than  in 
towns  to  the  north,  i.e.,  Danville,  Concord,  and  Antioch  (although  having 
a  focus  nearer  the  former).  At  the  more  northerly  area,  there  was  a  con¬ 
sensus  that  the  1/27  earthquake  was  felt  more  sharply  than  the  first  and 
there  was  confirmation  from  minor  building  damage  in  Antioch  in  the  1/27 
event  that  was  not  paralleled  during  the  1/24  earthquake  (Bolt  et  al. ,1981). 

Evidence  is  also  available  from  instrumental  readings.  The  first  set 
arises  in  terms  of  the  magnitudes  of  the  two  principal  shocks.  The  Wood- 
Anderson  seismograms  from  which  the  Berkeley  Richter  magnitudes  were  cal¬ 
culated  are  shown  in  Figures  4.11  and  4.12.  There  is  no  question  but  that 
the  maximum  amplitude  recorded  by  the  torsion  instrument  at  Berkeley  from 
the  1/27  earthquake  is  greater  than  the  maximum  trace  amplitude  from  the 
1/24  earthquake. 

At  the  Santa  Barbara  seismographic  station,  SBC,  to  the  south,  by  contrast 
the  1/24  earthquake  produced  substantially  larger  amplitudes  on  Wood-Anderson 
seismographs  than  did  the  1/27  shock  (see  Figure  4.13).  In  agreement  with 
the  oCher  evidence  for  strong  azimuthal  variation  in  radiation,  the  opposite 
ratio  occurs  at  the  northern  California  Areata  (ARC)  station.  It  will 
also  be  noticed  that  the  wave  duration  is  greater  in  the  case  of  the  first 
earthquake.  A  problem  for  comparisons  of  magnitude,  seismic  moment, 
intensity  and  so  on  for  this  pair  of  earthquakes  is  the  greater  complexity 
of  the  first  principal  shock.  This  is  evident  from  the  recording  of  three 
separate  peaks  (Figure  4.11)  indicating  multiple  events  superimposed  within 


Figure  4.11  Portion  of  the  Berkeley  100  *  magnification  torsion 
seismograms  showing  the  duration  and  amplitudes  in 
the  1/24  principal  shock.  Distance  between  time 
marks  is  1  min 
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Figure  4.12  Portion  of  the  Berkeley  100  x  magnification  torsion  seis¬ 
mograms  for  the  1/27  principal  Shock.  Same  scale  as  Figure 
4.11.  Notice  the  shorter  duration  and  larger  amplitude 
than  in  the  1/24  shock  as  recorded  at  BRK  about  30°  further 
west  of  the  NW  strike  of  the  ruptured  fault 


about  1  minute.  The  extended  duration  in  the  1/24  shock  has  two  conse¬ 
quences  of  importance  for  the  present  study.  First,  the  lengthier  shaking 
may  have  contributed  to  the  damage  reported  at  the  Lawrence  Livermore 
National  Laboratory  area  in  the  first  earthquake.  Second,  the  use  of 
duration  of  coda  length  for  discussions  of  directivity  focusing  can  be 
error-prone  in  the  case  of  multiple  shocks. 

The  mechanisms  of  the  ruptures  involved  in  the  two  principal  shocks 
can  now  be  summarized.  The  evidence  from  the  field  observations,  the  hypo- 
central  locations,  and  the  fault  plane  solutions  is  that  the  1/24  rupture 
began  to  the  north,  causing  right-lateral  displacement  along  the  northerly 
extension  of  the  Marsh  Creek  fault.  The  rupture  then  moved  along  the 
Greenville  strike  and  proceeded  for  about  20  km  southward,  producing  sur¬ 
face  ruptures,  with  the  faulting  ending  near  Highway  580  overpass  at 
Greenville  Road.  Such  unilateral  southern  progression  of  the  seismic 
wave  source  produced  radiation  conditions  which  could  have  enhanced  the 
strong  ground  motion  at  the  south  end  of  the  faulting  and  decreased  it 
somewhat  to  the  north.  The  position  of  the  focus  of  the  1/27  earthquake 
near  the  southern  end  of  the  faulting  suggests  further  southward  rupture 
of  the  Greenville  fault.  Field  evidence,  however,  indicates  no  surface 
faulting  south  of  Highway  580  along  the  Greenville  fault,  but  some  ad¬ 
ditional  displacements  along  the  surface  rupture  already  evident  after 
the  1/24  event.  A  plausible  explanation  is  that  the  second  earthquake 
involved  northward  rupture  velocity.  If  this  is  the  case,  this  event  is 
an  example  of  a  fault  dislocation  retracking  along  an  earlier  path,  either 
through  or  above  or  below  the  earlier  rupture.  The  depths  of  the  two 
principal  foci  of  11.8  and  14.5  km  suggest  deep  crustal  rupture  planes 
with  only  limited  expression  on  the  surface. 


The  radiation  of  ground  horizontal  peak  accelerations,  velocities 
and  energy  flux  in  the  1980  Livermore  Valley  earthquakes  has  been  subject 
to  a  close  analysis  in  a  recent  paper  by  Boatwright  and  Boore  (1982) . 

These  authors  examine  the  strong-motion  accelerograph  recordings  at 
27  stations  within  about  60  km  of  the  1/24  mainshock  and  the  1/27  shock 
for  systematic  variations  which  might  be  indicative  of  either  site  effects 
or  directivity.  They  find  that  the  variation  of  peak  accelerations  with 
azimuth  relative  to  the  source  is  reversed  for  the  two  events  and  agree 
with  the  intensity  information  that  the  mainshock  accelerations  are  larger 
to  the  southeast,  whereas  the  second  principal  shock  accelerations  are 
larger  to  the  northwest.  They  have  eliminated  the  site  effects  by  forming 
the  ratios  of  the  peak  accelerations  from  the  two  events  and  this  process 
indicates  that  the  source  directivity  in  this  case  caused  a  total  variation 
of  up  to  a  factor  of  about  5  in  the  peak  accelerations.  The  key  observa¬ 
tional  comparison  is  given  in  Figures  4.14  and  4.15.  They  find  that  the 
simple  application  of  the  directivity  function  given  by  equation  (3.1) 
fits  the  acceleration  data  well  for  a  reasonable  rupture  velocity.  They 
find  also  a  similar  but  smaller  azimuthal  variation  in  peak  velocities  and 
radiated  seismic  energy. 

It  should  be  pointed  out,  however,  that  the  effect  of  the  directivity 
in  this  case  is  being  measured  at  some  distance  from  the  source  itself  and 
with  accelerations  that  are  smaller  (generally  less  than  O.lg)  than  those  of 
engineering  interest.  In  the  1/24  shock,  all  but  one  station  had  a 
distance  greater  than  15  km  from  the  fault  source  trace  .  in  the  1/27 
principal  earthquake,  all  but  3  stations  had  distances  greater  than  10  km 
from  the  source  ;  in  this  case  two  additional  stations  recorded  at  distances 
of  4  and  8  km  from  the  source,  with  azimuths  of  246°  and  314  °,  and  peak 


Figure  4.14  Selected  SH  accelerograms  of  the  1/24  mainshock  superimposed  on 
a  map  of  the  Livermore  area,  showing  the  directivity  in  peak 
motions  and  relative  complexity  of  the  mainshock  accelerograms. 
The  acceleration  scales  have  been  adjusted  to  compensate  for  the 
expected  1/R  geometrical  spreading;  each  scale  has  100  cm/sec^ 
between  the  large  tick  marks.  (After  Boatwright  and  Boore,  1982) 
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Figure  4.15  Selected  SH  accelerograms  of  the  1/27  aftershock.  The  acceler 
ation  scales  have  been  adjusted  as  in  Figure  4.14.  Note  the 
strong  amplification  to  the  northwest,  relative  to  the  1/24 
accelerograms 


accelerations  of  0.26g  and  0.27g,  respectively.  Unfortunately,  the 
nearest  strong-motion  station  is  at  an  azimuth  which  would  have  high 
SH  motions  for  right-lateral  strike-slip  rupture  independently  of 
directivity  rupture  considerations.  The  recorded  data  also  show  a  few 
anomalous  points,  particularly  in  PGV.  The  data,  moreover,  do  not  directly 
measure  the  effective  directivity  in  a  narrow  lobe  along  the  fault  zone. 
They  do  show  that  a  moving  dislocation  effect  is  able  to  explain  the  ampli¬ 
tudes  of  predominately  SH  wave  motions  at  distances  greater  than  about 
25  km  better  than  a  static  (double-couple)  source.  The  evidence  is, 
however,  the  most  persuasive  to  date  for  significant  directivity  focusing 
contributions  at  distances  of  a  few  tens  of  kilometers  from  the  fault 
source  of  moderate  earthquakes. 


PART  V 


FINAL  ANALYSES 


5.1  General  Conclusions 

At  the  present  stage  of  strong-motion  studies  there  still  remains 
considerable  uncertainty  regarding  the  extent  and  effectiveness  of  direct¬ 
ivity  focusing  in  modifying  strong  ground  shaking.  Nevertheless,  a  number 
of  general  conclusions  can  be  drawn. 

First,  from  a  seismological  point  of  view,  the  effect  of  the  moving 
source  has  been  clearly  demonstrated  in  numerous  studies  using  seismo¬ 
graphs  located  at  both  moderate  and  great  distances  (i.e.,  the  far  field) 
from  the  source.  Such  studies,  however,  usually  concern  long-period  seismic 
waves  with  periods  above  2  to  5  seconds.  Secondly,  when  seismic  waves  in 
the  near  field  with  a  range  of  wave  frequencies  characteristic  of  engineered 
structures  (i.e.,  1  Hz  to  10  Hz)  are  considered,  there  is  as  yet  only 
limited  definitive  evidence  available,  and  this  is  somewhat  contradictory. 
Nevertheless,  the  likelihood  is  that  directivity  focusing  in  strong  motion  is 
limited  in  the  near  field,  but  that  the  elementary  predictions  of  the 
magnitude  of  the  effects  can  be  sometimes  modified  or  even  overwhelmed 
by  other  features  of  the  source  mechanism,  the  geological  variations  along 
the  wave  paths  and  within  the  fault  zone. 

Thirdly,  the  ratio  of  peak  horizontal  ground  motion  in  the  forward 
direction  of  fault  rupture  to  the  peak  motion  in  the  backwards  direction 
is  probably  greatest  for  ground  displacements  and  velocities  and  least 
for  peak  accelerations.  High-frequency  ground  accelerations  show  vari¬ 
ations  due  to  scattering,  attenuation,  and  source  asperities  that  mask 
directivity  effects. 


It  is  necessary  to  make  some  summary  comments  concerning  the  inferences 


of  Singh  (1981)  for  the  1979  Coyote  Lake  and  1979  Imperial  Valley  earth¬ 
quakes  (see  Section  4.1)  and  the  evaluation  of  Boatwright  and  Boore  (1982) 
for  the  1980  Livermore  earthquakes  (see  Section  4.3).  Singh  was  among  the 
first  to  demonstrate  from  strong-motion  accelerometer  records  the  reality 
of  directivity  focusing  in  the  geotechnical  context.  He  found  that  the 
effect  was  frequency- sensitive.  His  analysis  of  the  strong-motion  records 
available  to  him  suggested  that  peak  ground  acceleration  (high-frequency) 
was  affected  by  directivity  less  than  peak  velocity  (intermediate  frequency) 
and  peak  displacement  (low  frequency).  The  contours  of  peak  ground  accel¬ 
eration  for  the  Imperial  Valley  earthquake  were  complicated  and  did  not 
resemble  in  certain  ways  the  predicted  contours  for  a  simple  moving  source, 
and  he  pointed  out  that  the  peak  accelerations  obtained  at  two  stations 
closely  adjacent  to  the  Imperial  fault  were  opposite  to  that  predicted  by 
the  simple  model  of  directivity  focusing. 

By  comparison,  Boatwright  and  Boore  inferred  a  factor  of  5  in  accel¬ 
eration  focusing  in  the  azimuth  of  the  propagating  dislocation  along  the 
Greenville  fault.  It  has  been  pointed  out  (Section  4.3),  however,  that  in 
the  mainshock  crucial  strong-motion  instrumentation  was  not  available  close 
to  and  along  the  strike  of  the  Greenville  fault,  but  at  various  azimuths 
at  distances  beyond  15  km  from  the  seismic  source.  Boatwright  and  Boore 
were  therefore  required  to  interpolate  the  available  measurements  (see, 
e.g..  Figure  4.14).  Finally,  a  few  stations  in  the  case  of  the  Livermore 
earthquakes  did  not  have  recorded  accelerograms  in  agreement  with  the  over¬ 
all  radiation  patterns.  The  meaning  of  these  anomalies  from  the  predicted 
wavefield  has  not  as  yet  been  explored. 
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In  summary,  recordings  of  ground  motion  in  three  moderate  California 
earthquakes  of  1979  and  1980  provide  positive  evidence  for  a  measurable 
effect  of  dislocation  velocity  along  the  rupturing  fault  on  recorded 
strong  ground  motions.  The  consequences  are  superimposed  on  independent 
effects  arising  from  complex  source  mechanisms,  propagation  path  com¬ 
plications,  and  site  condition  variability. 


5.2  Engineering  Implications  for  Site  Evaluations 


As  summarized  in  Section  5.1,  the  evidence  presented  in  this  report 
indicates  that  directivity  focusing  of  seismic  waves  from  a  moving  fault 
rupture  has  been  confirmed  in  a  general  way.  The  main  effect  is  to  distort 
the  wave  radiation  patterns  that  occur  for  static  double-couple  force 
systems.  Because,  however,  most  strong-motion  parameters  are  averages 
of  samples  of  measurements  obtained  under  various  circumstances,  with 
variable  distances,  source  properties,  geological  conditions  and  azimuth 
from  the  seismic  source;  moving  source  effects  are  included  in  the  esti¬ 
mates  of  means  and  variances. 

Because  of  the  various  complexities  concerning  a  particular  site  and 
candidate  capable  faults  or  modes  of  faulting  in  its  vicinity,  quantitative 
incorporation  of  a  directivity  factor  in  site-specific  strong-motion 
parameter  estimates  cannot  be  recommended  at  this  time.  Although  there 
is  strong  evidence  from  the  field  instrumental  measurements  in  the  1979 
Imperial  Valley,  the  1979  Coyote  Lake,  and  the  1980  Greenville  (Livermore) 
earthquakes  for  directivity  focusing  effects,  the  actual  direction  of  the 
moving  source  cannot  be  predicted.  Therefore  a  site  may  be  on  the  high 
or  low  side  of  the  radiation  wave  field.  The  use  of  averages  of  all 
relevant  recorded  motions  will,  moreover,  always  incorporate  the  direct¬ 
ivity  effect,  at  least  in  the  sense  of  ground  motions  with  the  greatest 
expectation. 

For  example,  in  the  testimony  before  the  Atomic  Safety  and  Licensing 
Appeal  Board  on  the  Diablo  Canyon  Nuclear  Power  Plant  site  (see  Section  2.1) 
comments  were  made  on  directivity  focusing.  It  was  generally  agreed  that 
the  effects  of  focusing  are  generally  included  in  the  seismic  analysis. 


but  that  it  might  explain  some  observed  scatter  of  peak  acceleration, 
velocity,  and  so  on.  Furthermore,  the  concept  of  seismic  focusing  was 
characterized  as  generally  theoretical  in  substance  and  still  speculative. 
Testimony  indicated  that  the  effects  of  focusing  occur  in  every  earthquake 
to  some  degree  for  sites  where  there  is  propagation  in  the  direction  of 
the  instrument.  For  this  reason,  the  strong-motion  earthquake  recordings 
take  into  account  on  the  average  the  effects  of  seismic  focusing.  The 
difficulty  is  that  for  any  given  earthquake  on  the  Hosgri  fault  it  is  now 
not  possible  to  predict  whether  it  would  focus  energy  towards  or  away  from 
the  Diablo  Canyon  site.  There  was  testimony  that,  for  large  earthquakes, 
accelerations  may  theoretically  reach  2g  but  that,  to  date,  in  only  two 
cases  were  horizontal  accelerations  greater  than  lg  measured.  In  fact,  no 
witness  was  able  to  assign  a  probability  to  the  likelihood  that  given 
seismic  focusing  effects  would  be  associated  with  ground  motion  at  the 
Diablo  Canyon  site. 

Because  seismic  focusing  is  a  special  near-field  earthquake  ground 
motion  effect  which  may  be  experienced  at  a  site  towards  which  fracture 
propagation  progresses,  local  amplification  of  the  ground  motion  might 
arise.  The  effect  should,  therefore,  be  considered  in  shaking  estimates 
for  critical  structures  such  as  dams,  bridges,  power  plants,  etc.,  because 
such  focusing  can  significantly  increase  the  spectral  components  of  ground 
motion  employed  as  input  into  the  seismic  design. 

A  final  relevant  point  for  geotechnical  studies  is  that  modern 
assessments  of  near-field  strong  ground  motions  usually  include  elements 
which  generally  take  account  of  the  dynamics  of  fault  rupture.  It  is 
now  good  practice  to  include  at  an  appropriate  portion  of  a  near-source 


record  a  longer  period  pulse  which  corresponds  to  the  elastic  rebound  or 
"fling"  along  the  the  fault  as  the  dislocation  passes  by  the  site  (see 
Bolt,  1981).  This  fling  ensures  for  most  engineering  purposes  that  the 
longer  period  parts  of  the  response  spectrum  are  realistically  energetic. 


5.3  Future  Research 


The  general  survey  and  discussion  in  the  previous  sections  have 
brought  to  light  a  number  of  important  avenues  of  work  for  the  future. 
First,  intensity  studies  can  be  of  value  in  making  assessments  of  the 
practical  consequences  of  aspects  of  the  earthquake  generation  process. 

For  example,  the  detailed  field  evidence  of  intensities  along  the  ruptured 
fault  in  both  the  1906  earthquake  and  the  1952  Kern  County  earthquake 
along  the  White  Wolf  fault  are  most  important  correctives  to  the  view 
that  structures  are  most  heavily  shaken  along  the  fault  itself.  In 
future  earthquakes  where  surface  faulting  is  observed  and  the  requisite 
elements  for  intensity  assessments  exist,  every  effort  should  be  made  to 
make  detailed  accounting  of  damage  and  other  shaking  parameters.  In  this 
respect,  some  reassessment  of  the  calibration  of  the  Modified  Mercalli 
intensity  scale  is  probably  required.  In  particular,  it  would  be  most 
helpful  to  incorporate  more  of  what  is  now  known  about  the  effect  of  soils 
on  attenuation.  In  this  respect,  contributions  of  soils  engineers  in  re¬ 
assessments  of  intensities  and  reconstruction  of  the  intensity  scale  would 
be  valuable. 

Secondly,  it  is  clear  that  a  significant  advance  in  knowledge  of 
directivity  focusing  effects  depends  on  the  installation  and  maintenance 
of  strong-motion  accelerometers  along  active  faults  where  large  earthquakes 
might  be  generated.  As  the  studies  described  in  the  previous  sections  have 
shown,  there  are  still  only  a  handful  of  accelerograms  from  moderate  to 
large  earthquakes  that  allow  comparisons  at  the  front  and  back  of  a  moving 
source  dislocation. 

In  order  to  achieve  the  instrumental  requirements,  specially  designed 
digital  instrument  arrays  must  be  installed  along  active  fault  zones. 


These  arrays  can  be  designed  to  give  the  quantitative  detail  necessary 
to  model  properly  the  intricacies  of  strong  ground  motion  near  to  a 
source  moving  through  a  geologically  realistic  medium. 

On  the  theoretical  side,  the  discussion  has  demonstrated  that  more 
appropriate  theoretical  models  must  be  used  in  examining  the  quantitative  ef 
feet  of  directivity  focusing  in  the  frequency  ranges  of  interest  in  hazard 
assessment  and  for  engineering  design  purposes.  The  application  of  finite 
element  methods  to  the  complete  problem  promises  to  allow  the  introduction 
of  vertical  fault  zones  or  seams  in  the  country  rock  having  more  realistic 
seismic  properties.  In  the  meantime,  predictions  of  strong  ground  motion 
at  high  frequences  based  upon  simple  theoretical  models  that  do  not 
incorporate  adequate  variations  in  gradients  in  structural  and  elastic 
properties  and  wave  damping  must  be  viewed  with  caution. 
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APPENDIX  A 


GEOMETRICAL  CONSTRUCTION  FOR  RUPTURE  IN  A  LOW  VELOCITY  FAULT  ZONE 

When  a  vertical  (say)  fault  zone  is  taken  into  account  in  modeling 
the  seismic  wave  field  from  a  moving  dislocation,  computation  of  the  field 
becomes  more  complicated  than  in  the  simple  single  fault-plane  model. 
Solutions  for  the  seismic  wave  equation  for  boundary  conditions  appropriate 
for  such  a  geometry  do  not  appear  to  be  available  in  the  literature  (see 
Section  3.4),  although  some  finite  element  modeling  of  these  conditions 
has  been  started  by  the  author. 

It  is  of  some  value,  however,  to  give  an  elementary  picture  of  the 
effect  of  the  fault  zone  (or  seam)  by  means  of  a  family  of  wave  fronts. 
Consequently,  in  Figure  A.l,  a  modification  of  the  well-known  effect  of 
a  subsonic  unidirectional  moving  point  source  is  presented  for  comparison 
with  Figure  3.3.  (The  interesting  case  of  a  rupture  supersonic  in  the 
seam  but  subsonic  compared  with  the  wave  velocity  for  the  country  rock  is 
not  treated  here.) 

The  construction  in  Figure  A.l  is  purely  geometrical  and  does  not 
include  the  case  of  wave  conversion  at  the  boundaries  between  P  and  SV 
waves.  The  wave  front  can  be  taken  as  SH  waves,  for  example.  The  di¬ 
mensions  are  arbitrary,  but  the  wave  velocity  in  the  homogeneous  country 
rock  is  assumed  to  be  much  greater  than  that  in  the  homogeneous  seam. 
Refracted  waves  (D)  would,  by  Huygen's  principle,  radiate  in  general  from 
each  point  of  the  boundaries.  In  the  country  rock,  the  wave  fronts  (A1) 
would  move  ahead  of  the  corresponding  fronts  (A)  in  the  seam;  at  the  inter¬ 
section  of  the  front  with  the  boundary  there  would  be  a  conical  head  wave 
of  diffracted  type  (shown  as  a  dashed  line  D  for  the  initial  front  only) 


A1 


propagating  in  the  seam  as  a  leading  fringe  of  energy. 

A  few  inferences  on  the  likely  modifications  to  the  effect  of  direct¬ 
ivity  focusing  brought  about  by  the  seam  can  be  drawn  from  Figure  A.l: 

(a)  A  simple  quadrantal  pattern  with  lobes  would  no 
longer  describe  the  wave  amplitudes  in  all  the  near  field. 

(b)  The  concentration  of  energy  (focusing)  in  the  forward  direction 
is  not  as  great  just  outside  the  seam  as  inside  it.  Shaking 

in  the  seam  would  commence  before  the  arrival  of  the  original 
initial  pulse  producing  a  longer  duration  and  spreading  of 
the  wave  energy. 

(c)  Enhanced  damping  of  the  waves  in  the  seam  would  reduce  the  wave 
amplitudes  there  relative  to  waves  refracted  out  of  the  seam. 

(d)  The  diffraction,  multiple  reflection  and  progressive  interference 
effects  within  the  seam  would  be  a  function  of  the  effective  seam 
width  relative  to  the  wavelengths.  In  representative  geological 
conditions,  the  higher  frequency  (above  5  Hz)  waves  are  likely 

to  be  most  attenuated  by  the  multiple  interference  and  damping 
within  the  seam. 
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APPENDIX  B 


DIRECTIVITY  THEORY 


Consider  the  effect  of  a  finite  moving  source  on  the  seismic 
surface  wave  radiation  (Ben-Menahem,  1961).  The  solution  for  free 
surface  radial  Rayleigh  wave  displacements  for  a  vertical  strike-slip 
fault  is 


2  cos  0 


/  Kr 
8  o 


gr  (“) 


ei(0R-7r/4) 


(B.l 


where 

*R 


cob 


cos  0 

o 

-*R 


0q  =  azimuthal  angle  on  the  free  surface 
K  *  shear  wave  number  in  half  space 


r  =  radial  coordinate  on  the  free  surface 
o 


gr(“) 

b 


radial  layering  function  for  Rayleigh  waves 

dimension  of  fault  in  the  direction  of  motion 


a)  =  angular  frequency 

c  =  Rayleigh  wave  velocity 

K 

V  =  rupture  velocity 

r  ■  time 

sin  XR 

The  diffraction  factor  — - -  in  equation  B.l  represents  the  finiteness 
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of  the  source  in  the  direction  of  the  propagating  dislocation  and  its 
effect  is  to  produce  deviation  from  a  purely  dipole  radiation  pattern. 
This  function  produces  nodes  at 


Q 

•7^-  -  cos  9q)  *  mr  for  n  *  1,  2,  3,  ... 

R 

with  the  first  node  occurring  at  2  tt/u>  ■  b  [1/V  -  cos  0  / c  ] 

O  K 

Rewrite  equation  B.2  as 


—  c  -  cos  0  =  nX 

V  R  o 


where  A  =  2irV  /u>  . 

K 


This  condition  corresponds  to  destructive  interference  at  the  point 
of  observation  and  is  analogous  to  Fraunhofer  diffraction  through  a 
rectangular  slit.  Note  that  as  long  as 

nir  <  |  "Up  -  cos  0q)  |  <  (n+l)ir  (B.4 

R 

no  additional  zeroes  are  introduced  into  the  double-couple  pattern  and 
only  a  variation  in  amplitude  is  expected. 

Now  form  the  ratio  of  wave  amplitudes  at  equal  distances  from  the 
source  in  opposite  directions,  i.e.. 
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This  ratio  is  the  directivity  function  (Ben-Menahem,  1961).  The  function 
is  independent  of  the  source  time  function  and  is  dependent  only  on  the 
horizontal  dimension  of  source,  the  velocity  of  rupture  and  the  strike 
of  fault  with  respect  to  the  line  joining  the  station  with  the  epicenter. 
This  function  may  have  any  value  from  zero  to  infinity.  As  (b/X)  decreases 
it  approaches  unity  and  provides  a  measure  for  the  horizontal  dimension 
of  source  and  the  speed  of  rupture. 

For  a  simple  application,  equation  B.5  can  be  written  as 


(%  +  cos  eQ)  sin  X  "  COS  0°) 

D  "  7 - -  “ — ^b~T“ - —  (B*6) 

(-  -  cos  eo)  sin  -x-  (-  +  cos  0O) 

where  c  stands  for  either  Rayleigh  or  Love  phase  velocity. 

Both  the  directivity  function  (B.6)  and  the  phase  differential 
have  now  been  used  to  estimate  the  length  of  fault  rupture  and  the  rupture 
velocity  in  quite  a  number  of  earthquakes.  Source  parameters  are  obtained 
from  seismograms  from  two  stations  located  at  opposite  ends  of  the  fault 
or  from  surface  waves  of  different  order  (Rl,  R2,  ...  etc.)  recorded  at 
a  single  station.  If  the  waves  travelling  in  opposite  directions  do 
not  travel  the  same  distance,  the  seismograms  must  be  adjusted  to  equal 
seismic  paths. 

In  one  of  the  first  successful  regional  studies  for  California  earth¬ 
quakes,  spectral  nodes  for  Love  waves  from  the  1966  Parkfield  earthquake 
observed  at  Berkeley,  California  (at  a  distance  of  270  km)  were  used  by 
Filson  and  McEvilly  (1967)  to  infer  its  source  parameters.  Using  the 
first  node  at  T  ■  22.5  seconds,  they  estimated  a  lupture  velocity  of 
2.2  km/sec. 
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